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ABSTRACT 


Transient and steady state nucleate boiling in saturated LN,, and Fll 3 
at standard and near zero gravity conditions are investigated for the hori- 
zontal up, vertical and horizontal down orientations of the heating surface. 

A semi-transparent thin gold film on pyrex glass is used simultan- 
eously as a heater and resistance thermometer. Thermal recordings and 
high speed motion pictures of the nucleate boiling process are synchron- 
ized and recorded simultaneously. Test surface and bulk liquid tempera- 
tures are determined with an uncertainty of less than -1.68°F. 

The near zero gravity condition ( a/g^ .004) is obtained by the drop 
tower method. Heat flux is varied from approximately 300 to 30,000 
Btu/hr-ft 2 . 

Two distinct regimes of heat transfer mechanisms are observed during 
the interval from the step increase of power input to the onset of nuc- 
leate boiling; the conduction and convection dominated regimes. The time 
duration in each reqime is considerably shorter with LN^ than with F113, 
and decreases as heat flux increases, as gravity is reduced, and as the 
orientation is changed from horizontal up to horizontal down. 

In transient boiling, boiling initiates at a single point following 
the step increase in power, and then spreads over the surface. The de- 
lay time for the inception of boiling at the first site, and the velo- 
city of spread of boiling varies depending upon the heat flux, orienta- 
tion, body force, surface roughness and liquid properties, and are a con- 
sequence of changes .in boundary layer temperature levels associated with 
changes in natural convection. Following the step increase in power 
input, surfece temperature overshoot and undershoot occur before the 
steady state boiling temperature level is established. 
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The surface superheat required for the inception of boiling is 
found to be independent of heat flux and orientation. It is lower at 
a/gFO than at a/g a l and is influenced by surface roughness, being higher 
for smoother surfaces. 

The steady state nucleate boiling curve of heat flux versus surface 
temperature is found to be dependent on gravity and the orientation of 
the surface with respect to the gravity vector. The influence of gravity 
and orientation on the active site anl average bubble population densities 
appear to be quite small. 

The maximum bubble size and frequency of bubble departure are found 

to vary, depending upon the heat flux and orientation, according to the 
1 to - 

relation (D ) f = constant, 
max 

The natural convection and nucleate boiling components of heat flux 
are separated by using the concept of the "area of influence", and the 
heat transfer associated with each departing bubble is determined from 
the experimental measurements. The results are compared with theoretical 
predictions of the Latent Heat Transfer and the Vapor-Liquid Exchange 
models. 

This report was also a dissertation submitted by the first author 
in partial fulfillment of the requirements for the degree of Doctor of 
Philosophy in The University of Michigan, 1973. 

This report, along with Ref. [29], together constitute the total 
final report on work conducted under the subject contract. 
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CHAPTER I 


INTRODUCTION 


A. PURPOSE 

Nucleate boiling is one of the most efficient ways of transferring 
heat from a solid surface to a liquid, in that the temperature differ- 
ences are relatively small for large magnitudes of heat flux. Advancement 
in space technology has increased the importance of understanding nucle- 
ate boiling at zero gravity as well as at standard gravity. Incipient 
and nucleate boiling have significance in the long-term space storage of 
cryogenic liquids. The onset of nucleate boiling, called Incipient 
boiling, will determine the degree of superheat present in the liquid. 
Large liquid superheats constitute highly unstable systems and can re- 
sult in undesirable effects such as sudden and large rises in pressure 
when vaporization begins. Because of its importance, many aspects of 
nucleate boiling have been extensively studied in order to learn more 
about the basic mechanism Involved. Due to Its extreme complexity, 
however. It remains poorly understood, as is demonstrated by the diffi- 
culty of predicting nucleate pool boiling heat transfer adequately for 
various surface-liquid coafclnatlcns. A number of models have been pro- 
posed, but each has limitations. Many of the experimental results are 
somewhat Incomplete and conflicting. Many experiments involved the use 
of small wires, so that the test surfaces were small compared to bubble 
diameters, raising the question as to whether the results obtained are 
representative of behavior with large surfaces. 

The purpose of this work is to investigate the behavior of certain 
parameters in transient and stea<y-state nucleate boiling at standard 
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and zero gravity. Incipient boiling can be viewed as a transient process 
In Itself, since a change in the basic mechanism Is taking place. It was 
thus deemed pertinent that a study of the transient boiling process 
should be Included. In space applications in particular, where finite 
systems are present, all phase change processes encountered will be tran- 
sient. The onset of non-boiling natural convection may be anticipated 
to Influence the inception of boiling, as will the conduction processes 
in both the fluid and the heat transfer surface. A complete description 
of transient boiling must therefore Include these mechanisms. 

The fluids used in this study are liquid nitrogen and Freon 113 with 
three different orientations of the heat transfer surface relative to 
the prevailing gravity field; vertical , horizontal up, and horizontal 
down. The boiling point of Freon 113 is 117°F at atmospheric pressure, 
making it convenient to use at ambient conditions. Liquid nttrogen is a 
reasonably typical cryogenic liquid whose attribute of inertness makes 
it convenient to use. In addition, both fluids are electrical non- 
conductors, a necessity for the experimental technique used for heating 
and for surface temperature measurement. The three orientations of hori- 
zontal up, vertical, and horizontal down constitute the limits commonly 
encountered in practice. It is recognized that orientation has no sig- 
nificance under true zero gravity, but true zero gravity does not exist 
in practice. In space solar radiation pressure produces body forces of 
about a/g » 10 , and in drop tower experiments residual body forces of 
about a/g * 10 are present. 

The objectives of the present work are stated more explicitly as 


fol lows’- 
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1. To determine the influence of heat flux, heater surface orien- 
tation and gravity on the steady state and t ansient nucleate boil- 
ing heat transfer processes. 

2. To Investigate the role or significance of the bubble micro- 
layer In boiling heat transfer and bubble growth rates. 

3. To investigate the Influence of heat flux, heater surface or- 
ientation, and gravity upon the active nucleatlon site density, 
bubble population, frequency of bubble emissions, and maxi mum 
bubble departure size. 

4. To test some of the existing theoretical correlations with the 
experimental data obtained. 

B. LITERATURE SURVEY 

In this present work, both the transient and steady-state behavior 
of nucleate pool boiling at standard and near zero gravity condi tlons 
were studied. The effects of heat flux and orientation were also inves- 
tigated. 

The control and response of some thermal systems, particularly 
nuclear reactors, depend on the transient characteristics of the boiling 
process. Characteristics of boiling change considerably when very short 
transients are considered. It is therefore Important to know, for given 
fluids and heat transfer surface materials, the heater surface super- 
heat, delay time, and related heat flux needed to Initiate nucleate boil- 
ing. The literature on transient boiling heat transfer Is relatively 
limited. Those available are presented below, and then followed by a 
summary of the relevant points in each. 



Cole [1] and Rosenthal [2] Investigated transient heat transfer 
from nickel and platinum ribbons, respectively, to cubcooled water at 
atmospheric pressure. The test surface was used as a resistance thermo- 
meter to measure wall temperatures . Cole experimented with step power 
inputs while Rosenthal increased the power exponentially. A bank of 12 
volt D.C. batteries was used for the power source in both studies. 

Cooper and Lloyd [3] measured the surface temperature fluctuations of a 
heated glass plate in subcooled toluene at different pressures, as in- 
dividual vapor biibles form, grow and move off. The temperatures were 
measured by small thin metal and semiconductor thin films. Pressure was 
varied from 1-2 psia, heat flux was varied between 7500-15000 BTU/HR-ft 2 , 
and subcooling was varied between 0-15 degrees F. Garg and Patten [4] 
report on the temperature and pressure transients associated with bubble 
growth In and near the superheated layer around a 1/8“ diameter heated 
tube in saturated water. A 20-V, 400 amp, d.c. motor-generator unit 
was used to supply necessary power. The average surface temperature was 
deduced from measured stabilized bulk liquid temperature by a 35 SWG 
thermocouple. Hall and Harrison [5] Investigated the mechanism of 
boiling water with exponential 1y Increasing heat flux. A platinum ribbon 
was used both as heater and resistance thermometer. A set of storage 
batteries controlled by a transient power generator was used as the power 
source. Johnson et al [6] Investigated the transient pool boiling of 
subcooled water at atmospheric pressure. Energy was Increased exponen- 
tially in platinum and alimrtmsD ribbons, positioned both vertically and 
horizontally, and which were also used as resistance thermometers. The 
power source was similar to that used in Reference [5]. McCurdy [7] 
analytically studied the delay time for bubble formation In 
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waler moderated reactors, which were reported to be in agreement with the 
experimental result of Ref. [61. Lurie and Johnson [8] s fuelled transient 
pool boiling of subcooled water at atmospheric pressure, on a vertical 
surface with a step Increase In heat generation. An iron-nickel alloy 
ribbon was used as the heating surface and as a resistance thermometer. 
Power was supplied by a set of storage batteries. Klosterhouse [9] 
studied transient boiling using a small platinum wire as a heater and as 
a resistance thermometer. In saturated liquid nitrogen (lh 2 ) and liquid 
hydrogen (LH 2 ). A 12 volt D.C. battery was used for the step power 
input. Carpenter [10] made similar studies with a flat surface in 
liquid nitrogen, consisting of a pyrex glass substrate coate<' with a 
thin gold or nickel film and used as both a heater and a resistance ther- 
mometer. Graham [11] investigated the transient boiling of subcooled 
water and alcohol cm a horizontal chrome 1 ribbon cemented to a bakellte 
block. Surface temperatures were estimated from the readings of three 
thermocouples attached to the underside of the chrome! strip. A 12 volt 
battery was used to obtain the setp power increase for the experiments. 

With a step Increase In Imposed heat flux, what Is called a surface 
temperature overshoot was reported [1, 8. 11] , In which the heater sur- 
face tenperatjre Increases slcpilflcantly above the steady state tempera- 
ture corresponding to the particular heat flux, prior to the onset of 
nucleate boiling. Similar phenomena were observed with an exponentially 
increasing power input [2, 5, 6]. 

Once nucleate boiling is Initiated, generally at a single location. 
It spreads rapidly across the surface because of the presence of the 
superheated liquid. In so doing the liquid superheat Is reduced to near 
saturation condition, which can result In a momentary, decrease In the 



heater surface temperature below that of the steady tenperaturo corres- 
ponding to the heat flux. This has been observed with a ribbon and has 
been termed a temperature undershoot [11]. Owing to the rapid formation 
of vapor, more nucleatlon sites most likely have been activated than 
can be sustained In the steady state. This undershoot has alsn been ob- 
served in transient boiling with small wires [5], 

Time delays were observed between the Imposition of the heat fluxes 
and the inception of boiling [2, 6-11]. It was concluded [2, 8] that 
the transient process taking place before nucleate boiling began was 
solely a conduction process. As such it is possible to predict the sur- 
face temperature superheat by solving the equation of conduction. In 
cor-.arlng with predictions assuming a semi-infinite solid model, it was 
found that predictions were low In one case, [2], and high In another [8]. 
In another experimental work, using a Schlleren technique In conjunction 
with surface temperature measurements , it was observed that both trans- 
ient conduction and transient free convection processes occured prior 
to the onset of nucleate boiling, [11]. With a sufficiently high step 
increase In heat flux, the surface temperatures increased so rapidly 
that the process passed directly into and remained in the film boiling 
region [12, 13]. The surface temperature measurements of Ref. [4], while 
having a small response time, were highly localized and were Intended to 
show the effects of the growth and departure of the Individual bubbles. 

The present work Is Intended to experimentally determine the nature 
of the time delay, overshoot and/or undershoot In wall temperature, and 
the types of heat transfer processes taking place up to the establish- 
ment of steady state boiling, along with the effects of body force, or- 
ientation, heat transfer rates, and the surface characteristics on these 
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pnenomena. In addition, because of the experimental technique employed, 
it was possible to obtain data of Interest concerning certain aspects 
of steady nucleate boiling. 

Since it Is an efficient heat transfer process nucleate boiling has 
a wide range of application.*. Consequently a large nterfcer of research 
studies have been directed toward the different aspects of nucleate 
boiling. Many reviews of the general subject are available (e.g. Refs. 
[14, 15, 16]). Orly those works which have a direct relation to specific 
aspects of the present study will be reviewed here. 

Some existing experimental nucleate boiling heat transfer studies 
in reduced gravity fields are listed in Table 1, along with the signi- 
ficant parameters. \ number of works reported little or no influence of 
the reduced bod\ i orce on nucleate pool boiling process [17-23], On the 
other hand, others report a distinct decrease in the heater surface super- 
heat at a given heat flux as the body force Is reduced [24-28]. For 
the sake of convenience, the decrease in heater surface superheat for 
a given heat flux will be defined as a downward shift, and the opposite 
as an upward shift. 

The data of Clodfelter [27] Indicates a downward shift of 4°F. Papell 
and Faber [26] report about a 5°F downward shift. Littles [24] indicates 
an upward shift for vertical surfaces and downward shift for horizontal 
surfaces. Herte et al [25] observed a downward shift for horizontal up, 
an upward shift for horizontal down, and no change for vertical orien- 
tations, with reduced gravity. 

The seemingly contradictory results of Ref. [17-23] can be explained 
on the basis of the experimental procedures followed. Figure 10 in 
Ref. [19] Indicates an upward shift of 1.5°F for a/g - 0.20, but the 


Table 1. Nucleate Bolling Curve Shift from Standard Gravity to Near Zero Gravity 
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Authors report an uncertainty of 2°F in the data, rurthemore, the heater 
surface was spherical In shape, for which no preferred orientation of the 
acceleration vector exists. The Instrumentation in the experiments of 
Uslskln and Siegel [22, 23] was not sensitive enough to detect tempera- 
ture shifts less than 6 °F. There was a considerable scatter In the data 
of Sherley [18], most likely due to the lack of sensitivity In the tem- 
perature measuring technlqi-e. A least square statistical line fitted In 
the data of Ref. [18] by Littles [24] indicates a small shift such ,.<at 
the values of AT sat (Ty-T sat ) at zer0 9 rav ^ were °- 50|F smaller than 
the corresponding values at standard gravity. Indicating a slight down- 
ward shift. 

If nucleate boiling Is relatively Insensitive to buoyant forces, 
as reported by some of the investigators, and Is controlled by Inertial 
forces resulting from the dynamics of the bubble growth. It should also 
be Insensitive to the orientation of the heating surfaces even at stan- 
dard gravity. Marcus and Dropkln [30] noticed an Increase In heat flux 
for a given AT sat when they rotated the surface from horizontal up to 
the vertical position. Similar results were observed by Githlnji and 
Sabersky [ 31 ] when the orientation was changed from horizontal up to 
the vertical position, but they observed a considerable decrease In heat 
flux for a given AT when they changed the orientation of their thin 
ribbon surface from vertical to horizontal downward. On the other hand 
the highest heat flux for a given AT was reported with the horizontal 
down orientation In Freon 113 [24], Similar results were obtained In 
experiments with LN 2 and at standard gravity [25], 

Thtse latter results all Indicate that nucleate boiling i$ influ- 
enced by body forces. The magnitude of the effect of reducing body forces 
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must be investigated more thoroughly. A trend common In each of these is 
seen in Fig. 2 of Ref. [25], where at lower levels of heat flux the hori- 
zontal down orientation gives the highest heat flux for a given AT sat , 
and tite horizontal up gives the lowest, and the vertical lies in between. 

A possible explanation for this effec* lies in the role that microlayer 
evaporation plays in bubble growth. For the horizontal down orientation 
the buoyant forces tend to hold the vapor bi&bles closer to the surface, 
and thus could increase the effective area of the microlayer compared to 
the horizontal up and vertical orientations. Since the microlayer eva- 
poration produces large evaporation rates with small AT's, the heat flux 
for a given AT sat should be higher for horizontal down orientations. 

This s consistent with the observations of Hospeti and Mesler [32]. They 
report that the heat transfer associated with hemispherical bubbles is 
greater than that with spherical bubbles. Because of the influence of 
the buoyant force, horizontal down orientations tend to produce hemi- 
spherical bubbles while horizontal up orientations produce bubbles that 
are more spherical in shape. 

The behavior at higher levels of heat flux is opposite to that at 
It-, levels of heat flux in that the vertical and horizontal up orienta- 
tions give higher heat fluxes than the horizontal down orientations, for 
a given AT sat> as is shown In Fig. 2 of Ref. [25]. The explanation for 
this is probably related to an increase of the bubble population density. 
If the bubble population density is sufficiently high the interaction 
between the bubbles could reduce the relative effectiveness of the micro- 
layer mechanism. Oue to "crowding", the lateral bubble growth will be 
restricted, therefore the possibility of the formation of a microlayer 
could be reduced or eliminated. 
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To clarify these two pnenomena at the lower anu higher heat flux 
regions it is necessary to know how the active nucleation site density, 
maximum bubble departure sizes, and frequency of departures are influ- 
enced by gravity, orientation, heat flux, and the heater surface superheat. 

Many different models have been proposed to provide an explanation 
and a correlation for the nucleate boiling phenomenon. The apparatus used 
in the present investigation permitted simultaneous determination of 
heat flux, maximum bubble diameter, active site density and frequency of 
bubble departure for a given surface liquid combination. Therefore only 
the following two models which take full advantage of these results are 
evaluated here. 

1 . Latent heat transport model 

It was shown by Rohsenow and Clark [33] that for bubble radii great- 
er than l(f 7 ft the latent heat content, Qg, of a buoble of volume Vg 
is given, to a good approximation, by 

Q, = (l) 

with properties evaluated at the saturated state corresponding to the 
liquid pressure. Rallis and Jawurek [36] used Eq. {1} to write an ex- 
pression for the rate of energy carried away from such a surface by the 
bubbles, as 

^7 Va <« 

Equation (2) represents the latent heat transfer by the bubbles and 
strictly speaking applies only to saturated nucleate boiling. In terms 
of time averaged measurements and assuming bubbles are spherical. in 
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shape Eq. (2) may be expressed as 

<3) 

2. Vapor liquid exchange model 

Postulated by Forster and Grief <' 37 ] , it Is assimed that heat trans- 
fer in nucleate boiling takes place by virtue of the exchange of a volume 
of liquid at the bulk temperature (T^) for a volume of liquid at the wall 
temperature (T t ) each time a bubble departs or collapses. This model is 
expressed in terras of time averaged quantities by 

-7T = 0C f £rf**(TT-T h )-fcf 


( 4 > 


CHAPTER II 


EXPERIMENTAL APPARATUS 

The data obtained In the present study are heat flux, nucleation 
site densities, bubble growth rates and departure size, frequency of 
departures, and test surface superheats for the transient and steady- 
state nucleate pool boiling In saturated liquid nitrogen and in Freon- 
113 at standard and near zero gravity conditions. 

In order to study transient and steady-state boiling heat transfe* 
in cryogenic liquids at standard and zero gravity, the experimental appar- 
atus was designed and constructed to provide the following capabilities; 

1. Zero or near zero gravity conditions. 

2. A test surface large enough that It can represent actual 
surfaces, 

3. An optical and lighting system that permits continuous high 
speed photography of the heat transfer surface under all 
conditions. 

4. Instrumentation that permits recording the surface temperature, 
bulk liquid temperature , and heat flux rates accurately and 
continuously without any time delay for all circumstances. 

5. A system capable of synchronizing and sequencing 6-7 events 
In any desired order. 

6. A test dewar which can hold liquids having saturation tempera- 
tures ranging from -320*F to 120°F. 

7. A test package with little or no vibration as the body changes 
from standard to zero gravity condition. 


- 13 - 
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In the following paragraphs a detailed description of the apparatus 
used for the present study will be given. 

A. DROP TOWER AND DROP PACKAGE 

Zero gravity condition was obtained in the drop tower shown in Fig. 
1. It has a free fall distance of 3i ft (9.75 m) corresponding to a free 
fail time of about 1.34 sec. 

Figure 2 is a sectional view of the drop package and shows the de- 
celeration method used. The drop package had four different sections: 

1) release assembly, 2) piston assembly, 3) cylinder assembly, and 
4) the spike-metering pin assembly. The drop package al*~ serves as a 
container for the particular test vessel being used. It has space for a 
vessel or system having a maximum size of 19 Inches diameter by 36 Inches 
long. Prior to release the piston Is extended as shown in Fig. 2. At 
the end of the free fall period the cone and spike penetrate a barrel 
filled with sand to decelerate the cylinder. The piston continues to 
move, conpressing the air to a predetermined pressure corresponding to 
the maximum desired deceleration level. At this point a metering pin 
attached to the bottom of the piston asseably, which constitutes a vari- 
able area orifice, vents the air at a suitable rate to maintain the pres- 
sure approximately constant. At the appropriate position the metering 
pin closes completely to provide a final cushion of air. 

A Kistler model 808 piezoelectric crystal accelerometer was mounted 
on the drop package to monitor the deceleration levels. Figure 3 shows 
a typical result. The maximum deceleration was 25 g's, which was accept- 
able. 
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Figure 4 shows typical results of the acceleration or body force mea- 
surements (a/g) made directly on the test package during free fall. A 
Kistler model 303 servo accelerometer was mounted on the top cover plate, 
from which the entire assembly is initially supported. The vibrations of 
the cover plate-support assembly occuHng during release were picked up by 
the accelerometer. The mean readings of the test are plotted in Fig. 5 
and show a maximum value of a/g=0.QQ4 for the increase of the body force 
during free fall, due to air resistance. This is a 50% reduction of that 
reported in Ref. [29], and is attributed to the installation of an alum- 
inum fairing around the spike to give an aerodynamic shape. 

The elasped time for the free fall of 28 ft 7 inches was measured 
as 1.345 sec. The theoretical elapsed time of free fall in vacuum for 
the same distance is 1.335 sec. 

The calibrated zero shown in Fig. 5 arises from confuting the mean 
output signal between a/g=l and a/g= -1 by assuming the output from the 
servo accelerometer to be linear with body force, which is a severe assump- 
tion over such a wide range of output. The zero value measured by extra- 
polating to zero time is believed to be more realistic. 

B. TEST PACKAGE 

Figures 6 and 7 show the two different views of the test package. Fig- 
ure 8 shows a sectional view of the test package, which has six major 
components: 

1. Cover Plate: Because of weight limitations and the high vacuum 

sealing requirements for cryogenic temperatures, the 10 378 inch diameter 
inner section of the plate was made of stainless steel and the remaining 
of 20 1/2 inch diameter was made of aluminum. To reduce the vibration, 
the 1/2 inch thick cover plate is supported by two cross bars. The cover 



plate serves as a main supporting base for all of the major components of 
the test vessel, and mounts in the open end of the piston assembly in Fig. 2. 

2. Test Vessel Assembly: Two .065" thick wall concentric stain- 
less steel cylinders with 6.0" and 8.0" I.D.'s, closed at one end, each 
having four 2" dla, observation window, make up the test vessel assembly. 
These are bolted to the cover plate to complete the enclosure, and 1/2° 
thick pyrex glasses are used to make up four separate observation windows. 

In order to keep the heat transfer from surroundings at an acceptable 
minimum level when using liquid nitrogen, the space between the inner 

and outer cylinders is evacuated. The acceptable vacuum pressure for 
sealing was between 10~ 6 - 10" 4 mm hg. It was reduced to that level by a 
dual vacuum pump and continuously monitored except during the relatively 
short drop test. 

In order to seal the vacuum space at cryogenic temperatures , teflon 
coated stainless steel *'o a rings, .01" thick wall, were used at the win- 
dows and flanges. To compensate for any variations of the contractions 
between the parts when the dewar filled with Ut* two stainless steel spring 
washers were installed with each bolt. 

3. Test Piece Housing Tube: A 2-7/8" dia. flanged stainless steel 

tube with 0.065" thick wall was used to hold two relay mirrors and the 

test piece at the desired orientation. It is bolted to the cover plate 
and has the necessary openings for lighting and observation, test fluid 
fill line, and thermocouple and power lines. 

4. Optical System; A 200 watt mercury arc lamp, a spherical re- 

flectoi;, two lenses with 50 m focal length, and a mirror with 1/4 wave 

length flatness constitute the illumination side of the optical system. 
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In order to be able to illuminate the total or part of the test section 
at desired intensities, the system is designed such that the light source, 
each lens and the mirror can be moved independently. The precision 
movement necessary was accomplished by four sets of 1/4“ x28 N threaded 
rods and a sprocket chain assembly. Lenses and mirror were moved up and 
down inside a snugly fit 3-1/4“ dia., .065“ thick wall aluminum tube 
housing. The mercury lamp is located at the focal point of the 2" dia. 
spherical reflector, which Increased the available light intensity. 

The viewing side consists of a 1/4 wave-length flat mirror, a lens 
having 160 tun focal length and a camera relay mirror. In order to 
compensate for any misalignment during manufacturing of the optical system, 
the camera relay mirror has complete freedom of motion in every direc- 
tion. The viewing relay mirror and lens are made to move similarly to 
their illuminating side counterparts so that any desired magnification 
within the capacity of the system could be obtained. 

To Improve the stiffness of the total system and to provide that 
any vibrations occurring be in phase at the test surface, relay mirrors 
and lenses, the tubes housing the relay mirrors and lenses are bolted to 
the outer test dewar at the bottom. 

5, Camera Assembly: A Nova camera model 16-3 was used for taking 
high-speed motion pictures. It is designed to withstand up to 25 g 
ispact loads, and has a speed range from 100 pps to 10,000 pps In the 
16 tun format. The camera receives its power from a Model 1002 AC power 
supply. Camera speed is controlled by adjusting the voltage setting on 
the power supply. The camera speed/voltage setting was selected by 
reference to the speed charts provided by the manufacturers. Two neon 
timing lights within the camera were provided to record timing marks on 
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each edge of the film. One is operated at 60 cycles and its markings 
were used to calculate actual film speed. A 90 mm Mikrolor close-up 
lense model 16-508 is used as the camera objective lens. In some in- 
stances a 1/2" to 1" extension tube was used with the close-up lens to 
obtain increased magnification. 

The camera was mounted rigidly on a specially built support fix- 
ture, which was supported by two rails. The support fixture had the 
freedom to move in two directions in a horizontal plane along the rails 
and along the perpendicular direction to the rails. 

A Honeywell Pentax l°/20 exposure meter is used to determine the 
proper film exposure needed for different speed and lifting combina- 
tions. This meter was specially adapted for high-speed photographic 
works of up to 50,000 pps and has an angle of acceptance of 1°, which 
makes it ideal for the type of back lighting used in this work. 

6. Support Assent ly: The complete drop package is supported by an 
assent ly of three 3/16“ dia., 6x19 stainless steel ropes, a three-armed 
support fixture and a release rod. Stainless steel ropes connect the 
support fixture to the cover plate. Figure 6 shows the support assembly. 
The release assently has a circular groove In which balls are held in 
place by an external ring tapered on the inside. When the ring is lowered 
remotely by an air cylinder the balls rove out of the groove releasing 
the drop package, 

7. Drop Cable Assently: Three drop cables link the test package 
to the control area. The thermocouple outputs and test surface voltage 
connections are connected to the recorder via one of the drop cables made 
of eight 24 gauge copper and six 24 gauge constantan wires. This cable 
was double shielded to minimize noise pick-up. D.C. power to the test 

s 



-19- 


surface was supplied via a second cable. The third cable Is used to sudpIv 
the power for the camera, light source, and the synchronization signals to 
the neon timing lights. The cables may be seen In Figs. 6 and 7. 

C. TEST SURFACE 

Test surfaces used In studies of boiling heat transfer may be grouped 
Into two categories; 

1. Surfaces whose dimensions are in the same order of magnitude 
of the size of the bubbles growing on them. Small wires and narrow 
ribbons are examples of these types of surfaces. 

2. Surfaces whose dimensions are considerably larger than the 
bubbles growing on them. For example, discs, cylinders and spheres of 
relatively large dimensions. 

The first group of surfaces are generally used when surface temper- 
ature measurements are made using the surface as a resistance thermometer. 
High speed photographs of the bubbles mey be obtained simultaneously, 
since the small size of the heater surface usually will not generate 
enough vapor bubbles to obscure the field of view. The small size of 
this type of surface always raises the question as to how well It repre- 
sents the behavior of large size surfaces encountered In practical appli- 
cations. 

With the second group of test surfaces, the surface temperatures 
are generally measured with thermocouples somewhat removed spaclally from 
the actual heat transfer surface, giving rise to the necessity of extra- 
polating temperatures and Introducing time lags when dealing with trans- 
ient processes. When high speed photographs of boiling at the surface 
are desired, the field of view Is obscured by Intervening bifcbles. One 
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possibillty Is to use transparent surfaces such as glass or quartz and 
view the boiling process from underneath. This gives information as to 
the number of actives sites, the frequency of the bubble formation, and 
the departure sizes parallel to the plane of heating surface. A diffi- 
culty arises in transferring sufficient heat through the glass for boiling 
to occur because of the low thermal conductivity , and in measuring the 
temperature of the boiling surface. 

The first difficulty has been overcome In several research efforts 
[34, 35] by using a glass surface with a thin transparent layer of 
Sn0 2 baked on the surface which serves as the heat source when D.C. 
current is passed through it. This provides a large size surface with 
a uniformly generated heat flux at the surface. In one case the heating 
surface temperature was determined by measuring the lower surface tem- 
perature with an attached thermocouple and extrapolating to the upper 
boiling surface [34]. However, this technique will not give accurate 
results under transient conditions. 

In the present research, efforts were made to find suitable elec- 
trical resistive materials which could be deposited as a thin film on 
glass, which would be semi-transparent for optical viewing purposes, and 
which could serve simultaneously as a heating element and a resistance 
thei .nometer with negligible time lag. 

A number of materials, Including nickel, Sn0 2 , platinum, aluminum, 
and gold were Investigated. The most Important properties necessary to 
assessing the material which will serve as a resistance thermometer and 
heater are as follows: 



1. Temperature coefficient of resistance (TCR oroc ) defined by: 


with 


R(T) = R e f~ | +aC (T- To) 

- _ 



d R 
dT 


(5) 


where Ro ° resistance at the reference temperature T Q . 

2. Temperature coefficient of expansion (TC£ or /3 ) defined by. 

un = u [i +/9(t-t.)] (6) 


with 


P 


1 d L 
to d T 


where L Q Is the length at the reference temperature T Q . 

3. Transparency as a function of thickness. 

4. Adhesion of the film to glass. 

After extensive literature studies and testing, gold was selected 
as roost suitable. It is stable, free from oxidation, readily vacuum de- 
posited with reproducable resistance, has a relatively high temperature 
coefficient of resistance (TCR), and Is sufficiently transparent. 

Coming Pyrex glass surfaces were coated with a semi-transparent 
thin layer of gold and used as the test surface. The pyrex surfaces 
were first coated with &i approximately 300-400 A, thick gold film by 
vacuum deposition, and then the ends were coated with a relatively thick 
layer of silver as a contact area for the power supply and voltage taps. 
Figure 9 shows the overall dimensions used here. Thickness of the tlln 
film was monitored during the depositing process by means of a mass sen- 
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sltive crystal oscillator mounted adjacent to the test surface receiving 
the same rate of metal deposition as the glass surface. The measured 
natural frequency of the crystal was then a function of the metal thick- 
ness. This device was quite sensitive on a relative basis but did not 
provide an accurate measure on un absolute basis. Therefore several 
surfaces with different thicknesses of gold film were made and tested 
Initially. The most suitable one was selected by comparing their TCR 
and transparencies. Thereafter the film thicknesses were maintained 
approximately constant on the various surfaces used. 

Figure 10 shows the design of the current and potential connections 
which was determined to be the most effective by a process of develop- 
ment. 

During preliminary calibrations, to be described shortly, It was 
observed that the resistance of the gold surfaces varied with time but 
that the slope dR/dT remains constant. The resistance change was attri- 
buted to an aging effect. Moeller's work [38] dealing with the use of 
gold as a resistance thermometer indicates that Increasing the aging 
temperature (or annealing temperature) makes the resistance of the film 
much more stable, and specimens in that work were furnace soaked up to 
1080°K for 24 hours, A similar effect was reported in Refs. [39 ,40]. Be- 
cause of the use of pyrex as a substrate In the present application, such 
a high temperature was not possible. After experimenting. It was found 
that by tempering the thin films at 275°C for 15 minutes before calibra- 
tion, the variation In resistance with time was reduced considerably. 

Even after tempering at an elevated temperature to accelerate the 
aging process to a stable level, some shifting in resistivity with time 
took place at any given temperature leva'*, always Increasing but at a 
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much smaller rate, In such a manner that the slope dR/dT still remained 
constant. Several effects are believed to account for the Increase In re- 
sistivity of the thin gold film even after the accelerated aging process. It 
Is possible that the aging temperature used and the time of exposure for 
these gold films were not sufficient to provide the stable level. The 
usage and handling of these surfaces probably were causing a reduction 
of the thickness of the gold film which would account for the Increase 
In resistance. It was found necessary to use an Intermediate th. metal- 
lic element between the gold film and the pyrex during the process of 
fabrication In order to Improve the adhesion. Impurities In the gold an w 
change In the properties of the intermediate material could also account 
for these shifts. 

The continuous aging process caused no difficulty since Its effect 
was compensated by appropriate experimental techniques , which will be 
described later. 

A total of about 30 test surfaces were used during this Investiga- 
tion, all demonstrating the same behaH'v i . .hat their resistivity 
changed with time but the slope dR/dT remained constant. 

Figure 11 shows two complete calibrations of the test surface PCG-15 
In LN^ separated by an eight-month period of time. Also shown on the 
same figure are two sets of single point calibrations made Immediately 
preceding and following a boiling test. The boiling tests lasted on the 
order of four hours, and no shift occurred during this period of time. 

The two sets of calibration made eight months apart were almost parallel. 
Similar results were obtained with the other test surfaces used In this 
work and are shown In Appendix 0, along with the calibration of several 
other surfaces, obtained In the development of the gold film technique. 
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Slnce the variation In resistivity of the thin gold films was such 
that dR/dT remained constant, and because a complete calibration In 
requires approximately four hours a complete calibration associated with 
each test was impractical. Once the complete calibration of a test sur- 
face had been made, it was deemed sufficient to make only a single point 
calibration prior to and following each test run. 

Using the slight difference in slope of the two calibration curves 
in Fig. 11, it is possible to compute the error involved in the tempera- 
ture measurements had only one complete calibration curve been mad . 
followed by a single point calibration eight months later assuring that 
the resistivity gradient had not changed in that time interval. It can 
be shown that the fractional error in the measurement of a temperature 
by such a procedure will be given by 


AT er 2Ls ( dR /dT) , 
ATm ( dR /dl\ 


(?) 


where AT error a " T actual 


° ^actual ” ^single point calibration 

(dR/dT) ? a resistivity slope as measured at what is 
c taken as a single point calibration 


(dR/dT) i n resistivity slope as originally measured. 


For the two slopes shown in the Fig. 11 the error will be 
^Xtrrcr _ f Q.0 0VZ67 _ |1 _ 0.0/50$ 

at» nraara J 


This means that the error will be 1.5* of the difference In the tempera- 
ture between that to be measured and that of the calibration point. For 
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a temperature 20°F above the calibration point there would be an error 
of 0.3°F. 

During this investigation none of the tests conducted with the same 
test surface were more than four weeks apart. In addition, for experiments 
lasting more than four hours, three single point calibrations were made 
one immediately prior, one after and one at the midpoint of the testing. 

In the calibration curve resistances are presented in terms of 
ohms/square, since this quantity Is proportional to the thickness of the 
film and is not a function of the geometry of the particular heating 
surface used. Writing the total resistance as 

R = e -wr i»> 

R - total resistance 
P 3 specific resistivity -fl- cm 
L a resistor length 
W = resistor width 
t = resistor thickness 

The resistance in -tt/square Is given by PI t in the expression above. In 
thin films, the resistance Is no longer necessarily proportional to the 
bulk resistivity divided by the film thickness, so the resistance in 
ohms/square provides a more meaningful description. 

During calibration of test surfaces at different temperature levels, 
the D.C. current levels were varied between 5 to 95 cna to assess possible 
effects of current and self-heating. The heat fluxes corresponding to 
this range are less than 2 Btu/hr-ft , which Is negligibly small. Fig- 
ure 12 shows the resistance of gold film surface PC6-6 for various cur- 
rents. Even though data in this and similar subsequent tests Indicated 
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that it was not necessary to vary the current in the calibration tests, 
as a standard procedure the resistances of the test surface for four 
different current levels ranging from 5 ma to 95 ma were measured and 
then averaged for each calibre 'ji point. 

During calibration the electric;*! resistance of the test surfaces 
were measured using the circuit shown in Fig. 13, with the holder in- 
serted in various fluids to provide different known temperatures. Test 
surfaces were calibrated in the vicinity of two different temperature 
levels; the saturation ten$erature of FI 13 at atmospheric condi ti cm {about 
117°F), and the saturation temperature of liquid nitrogen (about -320°F). 
Calibration for the FI 13 temperature level was done in a constant tem- 
perature oil bath. To obtain the desired constant temperature a Colora 
Ultra- thermos tat filled with silicon oil was used. Temperature level 
varied from 25° to 100°C. Temperature constancy of the bath for this 
range was *0.005°C. Liquid temperature was measured with a Leeds and 
North rup platinum resistance thermometer (serial no. 1639500) with a 
sensitivity of i-OO^C, along with Leeds and Northrup Muller bridge with 
an accuracy of *0.02%. Resistance of the test surface was determined 
from the voltage drop across the test surface, Ey, and across the 
. 50274JJ. standard resistance, E s , see Fig. 13. BothEy andE $ were read 
on a K-3 potentiometer. 

Calibration of the test surfaces at the LN 2 temperature level was per- 
formed in the apparatus used In the work reported in Ref. [41], The test 
pieces were inserted in a stainless steel dewar filled with liquid nitrogen, 
and, the liquid saturation temperature was Increased with an inserted 
heater. A calibrated 12 inch Heise Bourdon tube gage was used to measure 
the vapor pressure of the LN 2 . The 25 psig range was divided into 0.1 psi 
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steps and had a calibration accuracy of *0.025 psi. Barometric pres- 
sures were measured with a mercury barometer. Pressure in the dewar 
was controlled by a mercury switch with an accuracy of *0.02 psi. LM 2 
temperature varied from 136 to 170°R and was measured by a copper-con- 
stant an thermocouple calibrated earlier. Readings of the pressure gage 
are used to determine the deviation of the thermocouple readings from 
MBS standards. 

D. POWER SUPPLY 

Figure 13 shows the power and calibration circuit diagrams, and 
consists of a 0.01 ohm standard resistance, a 0-16 ohm stepwise variable 
resistance and an adjustable d.c. power supply, consisting of nine 12 
volt, 72 amp. hour batteries which can be switched to provide 12, 24, 
or 36 volts. For test surface calibrations, the same power circuit was 
used except that the 0.01 ohm standard resistance was replaced by a 
0.50724 ohm standard resistance, and the 0-15 ohm variable resistance 
was substituted by a 0-2400 ohm variable resistance. 

Since the requirements for steady current is so stringent with the 
use of thin films, any variation in resistance wltn temperature of the 
variable resistors which nay be acceptable In normal applications could 
not be accommodated here. Therefore the 0-15 ohms variable resistor 
elements were constructed of chrome 1 ribbon and kept Immersed in approxi- 
mately ten gallons of transformer oil in order to maintain a constant 
temperature at the different levels of current. The 0.50724 ohm standard 
resistor was also kept In the oil bath. The 0.01 ohm resistor has a 
very large surface area but was Installed in an insulated box to minimize 

the effects of temperature drifts in the room. To minimize contact resis- 
tances, large capacity copper knife switches were used throughout. 
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E. INSTRUMENTATION 

To calculate the instantaneous heat flux values at different tempera- 
ture levels, the voltage drop (Ey) across the test surface and (E s ) 
across the .01 ohra standard resistance, giving the current, are recorded 
on a Sanborn recorder. The therraocou, le outputs, which give the bulk 
liquid temperature, and the synchronization and timing signals were also 
recorded cm other channels of the Sanborn recorder. 

Since the thermoelectric power of the copper-constantan thermocouples 
varied from 16 /^/k at liquid nitrogen temperatures to 39 / 4</ /K at 
Freon 113 temperatures , and a .05 ohm change in surface resistance cor- 
responds to a 1°K change in surface temperature, it was necessary to use 
highly sensitive instrumentation. Figure 14 stows a view of the instru- 
ments, power supply, and synchronizing system. 

F. THERMOCOUPLES 

All the thermocouples used during the experiments and test surface 
calibrations were made of 30 gauge copper-constantan. Prior to use they 
were calibrated at saturation temperature levels of FI 13 and LN ? . For 
the calibrations above room temperature, the Leeds and Northrup platinum 
resistance thermometer described previously was used. The deviation* 
from the standard tables furnished by Leeds and Northrup was less than 
6^>/ in the 1900 range. This corresponds to a deviation of less than 
0. 1°K in the temperature range of 290-325°K. The procedure used to 
calibrate the thermocouples at the saturation temperature level of 
was Identical with the procedure used to calibrate the test surface at 
at the same teog>erature level. The deviation from NBS thermoelectric 

Deviation « E $tan( j art j " pleasured* 



-29- 


tables publications (123 -A, R-188) was about at 5734^ range, 

which corresponds to a deviation of about .22°K in the temperature range 
of 60-80°K. In determining the liquid nitrogen temperature, this devia- 
tion curve is used. 

Three thermocouples were used to monitor the liquid temperature. 

THC I was located about 10 cms below, THC II was located at the same 
level, 1/2 cm away, and THC III was located at about 10 an above the 
test surface. The thermocouple outputs were connected to a recorder via 
a drop cable consisting of 24 gauge copper wires and 24 gauge constantan 
wires. To minimize the noise pick up, especially when operating at high 
sensitivities, all thermocouple wires were double-shielded. 

G. POTENTIOMETER 

For the calibration of test surfaces, temperature and the Sanborn 
recorder, a Leeds and North rup K3 potentiometer with a Rubicon Model 3550 
photoelectric galvanometer was used. The potentiometer has an accuracy 
of *0.01556 + 0.5^^ . Therefore measured temperature uncertainties 
were ±.1°K and *.04°K In liquid nitrogen and Freon 113, respectively. 

To reduce the effect of any building vibrations, the amplifier for the 
null detector system was suspended from the ceiling with soft springs in 
the supports. A Honeywell series 3100 galvanometer Indicated the un- 
balance in the circuit. 

H. SANBORN RECORDER 

A 7700 series 8 channel Sanborn recorder was used to record the 
thermocouple outputs, and voltages Ey and E s in Fig. 13 as a function of 
time. For recording thermocouple outputs and E the Model 8803 high 
gain d.c. amplifiers were used. The chart speed could be varied from 
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0. 25 mm/sec to 100 mm/sec in nine steps. In these experiments a chart 
speed of 10Q mm/sec was used. One second timing pulses were provided by 
a timing marker. The high gain amplifier had a sensitivity range from 
1000 ,^/div to 1 / /v '/div. The thermocouple outputs were recorded 

at a sensitivity of 5 or 10 /^/div., and V $ was recorded at a sensi- 
tivity of 50 /^v/div. The non-linearity of the recorder specified by 
the manufacturer was .1% of the full scale output, that is *.25^^ , 
±.5Q^v and * 2 . 5 ^v at sensitivities of 5 /•v'/dlv, 10 /W/div and 
50 /fv'/div., respectively. The recorded non-linearity contributes an 
uncertainty corresponding to less than *Q.Q3°K in temperature readings 
and about *,0025 amps in current readings. The preamplifier was capable 
of following from IQ? to 90% a step change in 5 milliseconds. With the 
sensitivity of 50 y^/div used, the system would follow 400 ,* v /sec., 
which was more than adequate for the processes encountered during this 
work. 

1. SYNCHRONIZING AND SEQUENCING SYSTEM 

Both light intensity and camera speed are adjusted before each run, 
to compensate the effects of the variation of active site density, spread 
velocity, bubble frequency on film quality. Camera speeds were varied 
between 1750 to 5000 pps during the present Investigation. A 100 foot 
roll of film will cover a time interval between 2.71 and 1.1 seconds at 
these camera speeds. As indicated before the maximum time available at 
the 0-G condition was about 1.4 seconds. In addition to the limited 
times available for the transient processes, steady state boiling at 
standard and zero gravity are successively recorded on the same roll of 
film. It was thus necessary that the themal and photographic record- 
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inqs of the various events and the controlled variables be carefully 
synchronized. A highly accurate and versatile synchronizing and se- 
quencing system was designed and constructed. 

Details of the sequencing and delay circuits for the various events 
such as starting the camera, releasing the test package for zero gravity 
operation, pulsing the synchronization signal to the camera and to the 
Sanborn recorder, and starting the boiling process are shown in Fig. 15. 

The complete circuitry consisted of five sub-circuits which are connected 
by a 3 PDT relay switch. The first circuit is the camera power supply 
and camera and event synchronizer. Its function Is to start the camera 
and any desired event with a specified delay. The camera is started first 
and puts the 3 PDT relay switch In operation after a set time delay. The first 
pole of the relay switch, which is normally open, operates the test pack- 
age release mechanism solenoid. The voltage across the relay solenoid is 
recorded on the Sanborn recorder. The second pole of the relay, which 
is normally closed, turns off one of the neon timing lights In the high 
speed camera. Since energizing the test package release solenoid and 
cutting off the neon timing light occur simultaneously, these two signals 
are used to synchronize the Sanborn recordings with the frames on the 
high speed film. The third sub-circuit operates this synchronizing neon 
aiming light at a frequency of 1000 cps via the normally closed pole of 
the 3 PDT relay switch. The fourth sub-circuit is another time delay 
switch. This circuit has a nine step time delay ranging from 77 ms to 
275 ms, and controls the test surface power circuit. The total surface 
power circuit consists of a 0.01A standard resistance, a variable re- 
sistance {0-15-fl. ), an adjustable D.C. power supply (6, 12, 24 volt) and 
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the test surface. The voltage droD across the test surface, Rj, and 
across the standard resistance, R^, are recorded on Sanbon. recorder 
charts simultaneously. R T indicates the test piece surface temperature, 
and R s indicates the current in the circuit. 
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TEST CONDITIONS 


A. FLUIDS 

Two different liquids, liquid nitrogen and Freon 113, were used for 
the present research investigation. Liquid nitrogen was selected as a 
representative cryogenic liquid since it is chemically Inert, readily 
available commercially, and has properties which are well known. 

Freon 113 was chosen as the test fluid for experiments in the vici- 
nity of ambient temperature levels. It is electrically non-conductive, 
reasonably chemically inert, considerably less hazardous than methanol 
or carbon tetrachloride, and its properties are reasonably well estab- 
lished. 

B. PRESSURE 

Most of the experiments were conducteu with the liquid close to or 
at the saturation state. Due to Interference by intervening bubbles in 
the optical path It was necessary to conduct several of the experiments 
with the liquid socv hat subcooled in order to obtain better photographs, 
particularly for steady conditions with the horizontal-up orientations. 
The saturation temperature at the heat transfer surface was obtained 
from available vapor pressure data for each fluid, using pressures cal- 
culated according to 

P =■ * Pvwp or (10) 

All tests were conducted at atmospheric pressure so that in general 
P °P 

vapor atmospheric. 
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C. LIQUID TEMPERATURES 

Bulk liquid temperatures were measured continuously by three copper- 
constant thermocouples at different vertical locations in the liquid. 

The temperature Indicated by the thermocouple located at the same level 
but removed horizontally 1/2 inch aw.y from the test surface was used as 
the bulk liquid temperature. Prior to each run with FI 13, two 200 watt 
immersed heaters were used to increase the liquid temperature to Its sat- 
uration level, nominally 117°F. 

D. HEAT FLUX 

For the steady state experiments four or five levels of heat flux 
were used at each orientation. Four different levels were selected so 
as to give a reasonable spread over a range, from the minimum level where 
boiling starts to a level close to burnout or to the heat flux where the 
bubbles are too numerous to permit adequate clarity of viewing. For 
some cases a fifth level of heat flux was selected such that boiling could 
not start, but could be sustained once Initiated at a higher level of 
heat flux. In other words, after boiling was initiated the heat flux 
was reduced to a level where boiling was just barely continuing. 

For the transient experiments three different levels of heat flux 
were used at each orientation. The lowest level was selected as that 
heat flux which would result in complete coverage of the test surface 
by boiling sites within 1.4 seconds of switching the power on. This 
time is dictated by the zero gravity period available in the drop tower. 

A second heat flux was selected close to burnout, and a third one at 
an intermediate point. 
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E. TEST SURFACES 

Three different test surfaces were used during this study because 
of unplanned burnouts on two of the surfaces. As a result several of 
tne tests were repeated with the different surfaces In order to observe 
If the nominally similar surfaces produced different results. This 
might occur due to differences In surface characteristics. In addition, 
for some orientations tests were conducted at heat flux levels other than 
tnose indicated above, depending upon the results of earlier runs. 

F. ACCELERATION 

Two levels of accelerations were used, a/g=l and a/g a 0. As indi- 
cated previously. It Is not possible to claim that true zero gravity 
exists in the tests conducted here. Based on the measurements shown In 
Fig. 5, It is possible to state only that the body force present on the 
drop package increases after release from some small level to a maximum 
value of a/g=.004 during the 1.4 seconds. 

Due to the short time duration at a/g=>0, nominally 1,4 seconds. It 
was not possible to determine decisively whether true steady state boil- 
ing existed at the end of the period, regardless of whether boiling 
existed at the beginning of the free fall period or was initiated during 
the free fall period. Two types of disturbances are used in this re- 
search which can result In transients: a step change in heat flux and 
a step decrease in body force. 

The transient arising from a reasonaoly large step change in heat 
flux is quite obvious, both at standard gravity and at reduced gravity, 
and these processes are referred to as transient boiling processe-. For 
the purpose of identification the process of a step reduction in body 
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force taking place with prior boiling Is referred to as steady state 
boiling at zero gravity, even though this process is also transient In 
the strict sense. This terminology is justified on the basis that a 
seemingly steady value of temperature was present within a short time 
after release of the test package. V.’ 'ether these results would corres- 
pond to those obtained with long term reduced gravity, say on tne order 
of hours duration, can be answered only with experiments of long dura- 
tion at reduced gravity. 

A summary of the ranges of the variables are given in Table 2 • 





CHAPTER IV 


TEST PROCEDURE 

Preliminary tests were conducted to establish a/g, the body force 
level on the test package during free fall, the deceleration level upon 
impact at the end of the free fall, and to check the optical and ther- 
mal recording systems. The tests for a/g and deceleration levels were 
performed only once, but the optical and thermal recording systems were 
checked before every run. 

A. FILLING THE DEWAR 

After preliminary tests were completed , the dewar was filled through 
specially provided lines. The depth of liquid in the test vessle was 
needed in order to calculate the hydrostatic pressure at the test sur- 
face which In turn was used to obtain the saturation temperature. 

The depth of liquid nitrogen above the test surface was measured 
approximately, by means of four resistor type liquid level indicating 
sensors located at different levels. The locations were such that it 
was possible to determine the LN^ level within 2.5 inches, and it was 
always kept between 15-17.5 inches above the test surface. 

The depth of Freon 113 was calculated by measuring the amount of 
Freon 113 put into the dewar before the start of a particular series 
of experiments. Depending upon the orientation of the test surface the 
dewar was filled with 4 or 5 liters of FI 13, whiwi gave between 5-7“ of 
liquid height above the test surface. At the end of the experiments the 
amount of liquid remaining in the dewar was measured. The variation in 
liquid level for any series of runs was always less than 1". 
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B. CALIBRATION 

Once the dewar was filled, a calibration procedure was followed 
prior, during, and following each test. This included calibration of the 
thermocouples and the test surface at the prevailing temperature of the 
liquid, and of the recorder. If the series of experiments took more 
then two hours, single point calibration of the test surface was per- 
formed three tiroes; prior to, at the middle and at the end of the experi- 
ments. For series lasting less than two hours, single point caliDration 
of the test surface was made twice; immediately before the start of and 
after the end of the experiments. 

The output of the thermocouples in the saturated liquid was measured 
with a potentiometer and used to find the deviation of the thermocouple 
output from the standard tables, assuming the deviation of the thermo- 
couple outputs to vary linearly from the liquid temperature to the refer- 
ence temperature (ice bath). 

For each recorder channel iced, the amplifier gain was calibrated 
by supplying a known input from the potentiometer and observing the cor- 
responding deflection on the recorder. After each amplifier was cali- 
brated, the appropriate zero suppression was applied again using the 
potentiometer as the signal source. After zero suppression the ampli- 
fier gain was checked once more. In addition, recorder outputs were 
checked against the potentiometer before and after each run. 

c - CAMERA AND LIGHTING ADJUSTMENT 

Prior to the start of a series of experiments, the test surface was 
illieninated at the desired Intensity by adjusting the light source and 
illiailnatlng lenses. Then the camera, relay mirrors, and viewing lenses 
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were adjusted for focusing and desired magnification. After focusing, 
the camera and relay mirrors were secured on their supports to prevent 
any motion due to impact and vibration. Then the liquid was heated to 
Its saturation temperature corresponding to the pressure at test surface, 
and boiling was Initiated at the des‘t.*ed heat flux level. While boiling 
continued, the light intensity at the test surface was measured with a 
light meter, and the correct camera speed was determined for an f/2 
aperture opening. The camera close-up lens was always set at Its maximum 
aperture opening ( f/2), vdiich provided the largest surface coverage for 
a given magnification. Then camera power supply voltage was set for 
the desired camera speed. Camera speed was not critical as long as it 
was greater than 1500 pps. 

D. SEQUENCING AND DELAY 

For a selected camera speed the approximate running time for 100 
feet of film was obtained from speed charts supplied by the manufacturers. 
For experiments at steady state conditions the event and camera synchron- 
izer, (see Fig. 15), was set to start the camera first* and then operate 
the test package release mechanism at about 1.4 seconds before the end 
of the film, thus making it possible to photograph both 1-G and 0-6 
boiling on the same roll of film. The timing signal generator was turned 
on, and neon lights in the camera were checked for proper operation. 

After settings, adjustments, and calibrations were completed, a 
trial run was conducted at a low recorder speed (5-10 ram/sec) to check 
for proper operation , including signal connections, sequencing and zero 
suppressions. The film was then installed in the camera, the observation 
window removed, the test package lifted and connected to the release 
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mechanism, the liquid brought to Its saturation temperature, and power 
turned on manually with enough time allowed to establish steady state 
boiling. Steady state test conditions were considered to have been 
reached when the outputs of Ey and E $ reached a constant level and 
remained there. At this point the one second timer was turned on, 
recorder speed was increased to 100 mm/sec, and the event synchronizer 
was put in operation with a remote control button, sequencing the events 
as planned. The camera starts first, and at an appropriate time the 
release mechanism is actuated, which releases the test package for 0-0 
condition with about a 10 millisecond delay. The Instant of release Is 
marked on the recorder chart by a sifpial provided by a dry cell In the 
package release circuit, and a corresponding signal is recorded on one 
of the film edges by turning off one neon signal light. 

For experiments with transient boiling, the synchronizer operated 
in a similar way with the additional function of turning the power on 
according to the pre-set sequencing schedule. The recorder is started 
manually, and the camera and event synchronizer Is activated by remote 
control , which In turn starts the camera, releases the test package, 
registers the synchronizing signals on both the film and reorder, and 
finally turns the power on at a set time. 

The release signal cm the recorder and the synchronization signal 
on the film are recorded about 10.3 milliseconds apart, and are used to 
synchronize the film frames with the thermal recordings. 

After the run was completed at impact the package was hoisted back 
into position, and the film removed and marked with an Identifying nw- 
bar. For each run the above operations were repeated. 


t* 



CHAPTER V 


DATA REDUCTION 

In these experiments, as in all measurements, errors are associated 
with determinations of test surface temperature, bulk and saturation liq- 
uid temperatures, heat flux, site density, bubble frequency, and the 
maximum bubble size. The detailed calculation of the related uncertain- 
ties in these quantities is given in Appendix B. In the following sec- 
tions the methods used for calculating these quantities are described. 

A. TEST SURFACE TEMPERATURE (Ty) 

Test surface temperature was determined from the measured test sur- 
face resistance, Ry, the single point calibration resistance R c » and the 
corresponding temperature T c , and the slope of a calibration curve, 
®=dT/dR, according to the equation 

T t - Tc i tn (R-r - R c ) <11) 

The values of Ey and E & were obtained from the recorder charts using 
appropriate zero suppression and the sensitivity value at which the re- 
corder was operated. 

The level of the zero-suppression applied was varied for each test 
dependent upon the Ey and E g range and the sensitivity of the recorder. 
Sensitivity of the recorder was selected such that all variations during 
the various types of experiments (transient, zero gravity, etc.), could 
be recorded in the linear range available on the recorder chart with a 
minima of noise pickup. Test surface voltage drops Ey were recorded by 
a low gain an^llfier, at a sensitivity of 0.005 volts/div. and standard 
resistance voltage drops E £ were recorded by a high gain amplifier, 
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wlth 50 yj.v/di v. sensitivity. The standard resistance R $ was .01 ohm 
for all cases. 

As explained In Appendix B, the values of Rj. , R c , T c and dT/dR were 
determined with uncertainties of *.051%, *.033%, *.028% and *1.46%, re- 
spectively. Therefore, wall temperatures were determined with an uncer- 
tainty of *1.33°F for the worst case with F113. For liquid nitrogen the 
uncertainty was less than *1.68°F for the worst case. 

It should be emphasized here that since the resistance of the total 
surface area is measured, the corresponding temperature represents only 
the mean temperature over the entire area. This mean will differ from the 
representative surface temperatures only in those cases where highly 
localized differences in behavior occur, such as during the spreading of 
nucleate boiling across the surface. 

B. BULK LIQUID TEMPERATURE < T /> 

Bulk liquid temperature was measured by copper-cons tan tan thermo- 
couples, using both the £3 potentiometer and the recorder. The Sanborn 
recorder was operated at a sensitivity of 5 yV-v/div. The maxinura uncer- 
tainty in the bulk liquid temperature was *.076% corresponding to less 
than *.09°F for the worst case. 

C. SATURATION TEMPERATURE IN THE TEST SURFACE (T $at ) 

The saturation temperature of liquid at the test surface was deduced 
from the pressure at that location. Pressure at the test n-rface was 
obtained by adding hydrostatic pressure at the test surface to local at- 
mospheric pressure. The main source of error in determining saturation 
temperature came from uncertainty in the liquid depth. The uncertainty 
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In the liquid level was less than *1 Inch for F113 and *2.5 inches for 
liquid nitrogen, corresponding to *.1S°F for FI 1 3 and *.05°F for L^. 


0. HEAT FLUX 


The heat flux to the liquid at the heat transfer surface in the 
steady case was calculated from the total power generated on the surface, 
subtracting the heat transfer to the glass substrate. 


(q/a\ = (a/4) T -(a/Ak„ 


(12) 


Since the test surfaces were kept in a snugly fitted teflon holder, heat 
losses through the edges are neglected. For steady state runs heat trans- 
fer to the glass substrate In FI 13 and Lf^ is obtained from the curves 
given In Fig. 16 and 17, respectively. These curves were generated from 
the solution of the heat transfer problem shown on than, involving 
boiling, conduction, and convection. The uncertainty associated with 
total heat flux measurements were less th*n *1Z. The uncertainty in 
values of (Q/A) -| oss » comes from the imcertainty associated with determina- 
tion of wall temperature and uncertainty associated with natural convec- 
tion correlations, the latter contributing the most. 


E. ACTIVE SITE AND POPULATION DENSITY 

The active site density and average population density were deter- 
mined from the 16 on high speed motion picture films, using a Vanguard 
motion analyzer. The active site density is the menber of active nucleat- 
ing sites per unit area present under the given conditions. The average 
population density is the average number of bubbles present on the sur- 
face per unit area at any given time. During the growth period, which 
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is also called the active period, a bubble will be present at the site. 

No bubble is present at the site during the waiting period, and thus the 
population density is always less than the active site density. After 
the film was threaded on the analyzer the synchronization signal was 
brought onto the screen and the frame counter was set to zero. The re- 
lease of the test package takes place about 10-11 ms after the release 
signal, corresponding to 20-55 frames for film speeds of between 2000- 
5000 pps. Pictures after these 2Q-55th frames pertained to zero gravity 
boiling, while all earlier ones pertained to standard gravity boiling. 

Each frame of the film corresponding to a/g°l was studied carefully 
by running tne film on the analyzer backward and forward until all the 
active sites for standard gravity were located. After site counting at 
standard gravity was completed, the zero gravity part was analyzed. Since 
bubbles stay on the surface and continue to grow under zero gravity con- 
ditions, they coalesce as they grow, covering some of the sites on the 
surface, and the short time at zero gravity makes it impossible to deter- 
mine all the active sites at zero gravity. Therefore only the newly ac- 
tivated sites were located and counted. The nunber of additional new 
sites activated during 0-G was added to the total of active sites counted 
during the 1-G condition. A new total was thus obtained ar^ used to cal- 
culate the active site density for the G-G conditions. Therefore, the 
active site density values for G-G are always larger than for 1-G con- 
dition. It is not possible to state that all the sites active at 1-G 
would continue to be active at 0-G, but as may be seen in the experimental 
results in Fig. 79, for example, the surface temperatures are higher at 
0-G than 1-G for a given heat flux, and it thus seems reasonable to assume 
that all the sites active at 1-G would therefore remain active at 0-G. 
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More than 502 of the surface area was photographed and enlarged 
more than 50 fold on the analyzer, and every frame of the 100 foot roll 
of film was viewed Independently. For one test run the film was analyzed 
by three Individuals at different times. The maximum variation In the 
active site density between the thret results was less than 72. It Is 
assumed that the portion of the surface photographed was representative 
of the total test surface. 

Population densities were obtained from 20 arbitrarily selected 
frames of the film. In each of these randomly selected frames the number 
of bubbles attached to sites were counted. The average of the randomly 
selected 20 frames was used to calculate the average population density. 
To find the representativeness of the sample size selected, average popu- 
lation densities were determined for 5, 10, 20, 30 arbitrarily selected 
frames from the same runs. The variation between using 5 and 10 frames 
was *6%, the variation between using 10 and 20 frames was *2.32, the 
variation between using 20 and 30 frames was *.752. In order to keep 
data reduction time as low as possible and still obtain accurate results, 
20 frames Is chosen as a basis to determine the population site density. 

F. BUBBLE FREQUENCY 

The period of the cycles between bubble emissions was determined 
by a detailed examination of the films projected on the analyzer screen. 
Depending upon the heat flux levels, the fraction of the total active 
sites that were used for determination of the bubble cycle period varied 
between 10-952. This variation was dictated by two factors. At high 
heat flux levels so many active sites existed that It became Impractical 
to follow each site. Secondly, the Increase in the number of total sites 
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with high heat flux reduces considerably the number of sites where the 
bubble growth period can be followed for at least one complete cycle, 
because of the periodic obscuring of portions of the surface. 

In the determination of bubble frequencies, the growth of a single 
vapor bubble at each site was followed frame by frame. The time inter- 
vals between bubble Initiation, departure, and the initiation of the next 
bubble at the same site were obtained. The elapsed time between initia- 
tion and departure of the vapor bubble was called the active period, and 
the time from departure to Initiation of the next bubble was called the 
inactive period. For each site analyzed, bubble growths are followed 
for at least 2 or 3 cycles until variations In the periods were reason- 
ably well determined. The average cycle length was used as the period 
for that site. The frequency of vapor bubble emission was obtained by 
taking the reciprocal of this quantity. Figure 18 shows a typical bubble 
appearance, growth and departure in LNg at standard gravity. The small- 
est diameter plotted corresponds to the frame in which the bubble first 
appeared. This means that the growth initiated some time in the inter- 
val of the previous frame. An inherent uncertainty thus exists in the 
origin of the growth period. This Is discussed In Ref. [46], 

G, BUBBLE SIZE 

The departure and the maxlrura diameter attained by the vapor bubble 
while attached to the heat transfer surface was measured , a tong with the 
bubble periods, from the films projected on the analyzer screen. The 
photographed area was enlarged 50-60 times on the screen, corresponding 
to a 7-8 power linear magnification. The film was again analyzed frame 
by frame and when the diameter of the vapor bubble attached to the sur- 
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face reached a maximum size It was measured In two mutually perpendicular 
directions by means of crosshairs on the analyzer, using a scale grad- 
uated In thousandths of an Inch. The average value of these two measure- 
ments was taken to be the maximum bubble diameter. Measurements of 
bubble sizes were taken for at least ? or 3 consecutive periods to deter- 
mine the variation In measured diameter, and to obtain better accuracy by 
reducing the chance of missing the maximum bi&ble size If It had been 
attained In the Interval of time separating two consecutive frames. 

Two 46 gauge wires, spaced nominally .1 Inch apart, were visible 
In the field o x view for every run. The precise spacing was measured by 
a traveling microscope each time the test surface was Installed. 

H. TRANSIENT PERIODS 

Two built-in neon lights were used to record timing and synchroni- 
zation marks on the edges of the film. In the transient runs, power was 
turned on automatically about 54.6 ms after the synchronization signal . 
was recorded, as Indicated by turning off one of the neon lights. The 
actual film speed was calculated from the recorded 60 cps signals on the 
second edge of the film. 

Figure 19 shows a typical representation of the transient surface 
temperature behavior occurlng with a step change In heat flux. The var- 
ious domains are believed to represent distinct physical phenomena taking 
place and are defined as follows: 

Tj - Is the time delay between the beginning of power 
input and onset of natural convection. It is thus the 
elapsed time In the conduction dominated regime. 
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^IB ' 1s the elapsed tine between the beginning of power 
Input and the Inception of the first boiling site as ob- 
served on the photographs. 

Tj - Is the elapsed tine In the convection dominated regime, 
and Is terminated when the surface temperature begins to de- 
crease because of the spreading of nucleate boiling. 

T, + T 2 - Is the elapsed time between the beginning of power 
Input to the maxi tarn surface temperature. 

T 3 - elapsed time from to the end of spreading of 
nucleate boiling over the surface. 

T s - represents the elapsed time from beginning to the end 
of the spreading of nucleate boiling over the surface. 

Ty * 7*2 + T 3 - represents the total elapsed time from the 
step power Input to the completion of the boiling spread over 
the test surface. 

T-j values were obtained from thermal recordings of Vy or from Ty 
vs time plots as In Fig. 20, which presents data obtained with a verti- 
cal orientation In F113 at standard gravity, for various levels of heat 
flux. Ti Is *n approximate value since the rate of increase In the 
surface Temperature Ty does not change sharply with the onset of convec- 
tion. It tapers down slowly. 

^IB values were obtained from the films. Since only 50-80* of the 
area of the test surface was photographed T IB values could not be ob- 
tained with the same accuracy for all cases. The syabol U IS“ in Fig. 20 
Indicates the point at which the first nucleating site is activated. 

When boiling was initiated outside the observed viewing area exact values 
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of 7 IB were not obtainable. Values for those cases are taken as the 
time of appearance of the first bt&ble, thus giving values of 7 ^ 
larger than may actually exist, 

T 2 values were not determined alone. : -stead the sibb of 7 ^ * 7 2 
values are measured. If 7 \ values were needed these can be obtained by 
subtracting 7 ^ from 7] * 7 * 2 * 7 ] + 7% were determined from the 

Sanborn recordings. Power was turned on 64.9 msec, after the Sanborn 
synchronization signal. Recorder speed was always kept at 100 em/sec. 

7\ + 7i values were obtained by counting the nunber of divisions from 
power on to the point that surface temperature reached a raaxlmun value. 

T 3 values were determined from the thermal recordings. It cor- 
responds to time between the time when surface temperature reaches 
its maximum value and the end of spread of boiling . T s values were 
determined from photcgraphi liras by observing the tin? between the 
beginning and the end of the spread 0 .' nucleate boiling. In addition to 
the spread times, T 3 and T s . the rate of boiling spread over the 
surface can be represented by the percent area covered by active sites 
versus time values and by deterralnlng the spread velocity. The fraction 
of the surface area covered by nucleating sites Is obtained by dividing 
the area covered by active sites at a given time by the total photographed 
area. Typical photographs of the spread of the active sites are shown 
in Figs. 95 and 96 of Appendix C. 

Spread velocity values, V s , are obtained in two steps. The dis- 
tance traveled by the bubble front from the starting point In two differ- 
ent directions were obtained and plotted against tloe. The average 
spread velocities are then obtained by averaging the slopes of these 


curves. 
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Total delay time values, T-\* * Tz * 7V " er - obtalr.ed by court- 
ing the timing signals from the frame where power was turned on to the 
frame whose surface was covered completely by the active bubble sites. 
This quantity was also Independently obtained from the thermal recordings 
by observing the variation of test surface terperature, as may be seen 
In Fig. 20. The variation between these two results for the majority of 
cases was not more than a few percent, with the temperature results 
giving larger values. This was reasonable since It was not possible to 
photograph the total area of the test surface, and the spread of boiling 
may have continued beyond the observed area. This unobserved part of 
the boiling was reflected in the thermal recordings. 



CHAPTER VI 


RESULTS AND DISCUSSION 

In this section, the results of the experimental measurements are 
presented and discu'sed under two main headings; Transient and Steady 
State Boiling Heat Transfer. These are presented in the form of graphs 
showing the influence of heat flux, orientation gravity and in some 
cases surface roughness on boiling heat transfer in FI 1 3 and LN~. An 
index of all tests conducted is presented as Appendix F. 

A. TRANSIENT BOILING HEAT TRANSFER RESULTS 

During the transient part of the process between the step power 
input and the point where the surface is completely covered by active 
boiling sites, three distinct regimes of heat transfer mechanisms were 
observed; the conduction, convection and boiling dominated regimes. 

The elapsed time in each regime varied depending on the heat fli'x, orien 
tation, body force and liquid used, which constitute the independent 
variables. The dependent variables are the transient test surface 
tenpe^ature, the various defined delay periods, the 9pread velocity, and 
the fractional coverage of the surface by nucleating sites. 

At the outset of this study, it was anticipated that the timewise 
activation of nucleation sites as the boiling spread over the surface 
would be determined. However, in viewing the films it appears that once 
nucleation takes place in a particular region, the boiling spreads by 
a process in which a vapor lawyer grows within the superheated liquid 
layer adjacent to the heating surface. Once the vapor front has passed 
a particular region, removing the superheat in the liquid, surface ten- 
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si on forces then result in the formation of b'/bbles and the normal 
nucleate boiling process follows. During the spreading of the vapor 
la}®r it was thus not possible to discern active nucleating sites. 

1. Test Surface Temperature Results 

Transients of the test surface temperatures are presented in the 
form of plots of surface temperature, Ty, versus time with levels of 
heat flux, orientation, body force, and the liquid used as parameters. 
Figures 20-23 and 24-27 depict the transient surface temperature results 
obtained with FIT 3 and LN^, respectively. The effects of heat flux, 
orientation, reduction of the gravity and different surfaces are pre- 
sented In these figures. Figure 28 shows the transient wall superheat 
for both FI 1 3 and LN 2 for the vertical surface at a/g*l. 

These all show a definite wall tenperature overshoot, AT q , as has 
been observed elsewhere [5,6,8-11]. The results reported in Refs. [8-11] 
should be considered most comparable since step Increases in power in- 
puts were used. 

It appears that the amount of the step Increase In heat flux has 
some effect on the degree of overshoot. Variation of AT Q with heat 
flux for vertical surface in FI 1 3 and LN 2 can be seen in Figs. 20 and 
24, respectively. These results along with others are replotted as AT Q 
vs. (Q/A)y in Fig. 29 for convenience, where the increase of AT q as 
heat flux Increases is more clearly noted. Similar results, thp increase 
of AT q with Increasing heat flux, have been observed In transient boil- 
ing studies with water [6]. 

On the other hand. Ref. [$] reports diminishing overshoot with 1n- 
crea: 4 ng heat generation rate. Also, Refs. [8] and [11] report no stir- 
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face temperature overshoot with saturated boiling, which Is contrary to 
the results of the present study and Refs. [9], [10]. These discrepan- 
cies could be due to the differences in surface roughness and heat capa- 
city of the surfaces used, and the sensitivity of the temperature mea- 
surements. As will be discussed iaUr, surface roughness has a pro- 
nounced effect on AT q : smoother surfaces result in a higher degree of 

overshoot In order for btijble nucleatlon to take place. Since no in- 
formation about the surface roucfinesses is given in Refs. [8,11], it is 
not possible to assess this effect. Furthermore, an insensitivity of 
about *5°K in temperature measurements was reported in Ref. [8]. A 
narrow chrrrael ribbon 1/16 inches wide cemented to a bakelite block 
served as the heater in Ref. [11], and temperatures were measured by 
three thermocouples attached to the underside of the chrorael strip. No 
information is given about the method of attachment and the accuracy of 
the temperature measurements , which give only local temperatures at the 
points of attachments. In addition to the possible Inaccuracies and 
time delays in the teaperature measurements, it was not possible to st&- 
stantiate a conclusion of “’no overshoot at saturated boiling 1 from 
Fig. 3 of [11], where transient surface temperatures for 45°, 23° and 6° 
sti>cooling are shown. 

As reported in Ref. [8], the lack of observation of any tenperature 
overshoot in Refs. [8] and [11] could be attributed to the heat capacity 
of the test surfaces used. Metal lie ribbons were used as a heating 
surface, with a volumetric heat generation cohered to a virtual plane 
heat generation with the thin gold film In the present Investigation. 
Thus, for a given power Input, the rate of temperature rise of the 
ribbons will be lower than that of the thin gold film even with the 
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presence of the glass substrate, which has a relatively low thermal 
dlffusivity. Transient effects in the liquid may thus be masked by 
transients in the heater surface, in the case of solid metallic surfaces. 

Transient surface temperatures with the vertical, horiaontal-up 

and horizontal -down orientations at a/g=*l are presented in Figs. 20, 22 

and 23 for F113 and in Figs. 24, 26 and 27 for IN, . Since different test 

surfaces are used in the experiments with FI 13, comparison of Figs. 20, 

22 and 23 to determine the effect of orientation on AT must be made 

o 

with caution. On the other hand, surface PCG 29 was used in the ex- 
periments with LN 2 for all orientations, and tne results presented in 
Figs. 24, 26 and 27 can be cos*) a red to determine the effect of orienta- ' 
tion on AT Q . For convenience, the AT Q 's are replotSed vs (Q/A)^ in 
Fig. 29. Figure 29 indicates that the effect of heat flux on AT Q is 
more pronounced with horizontal-up orientations and the vertical than 
the horizontal -down orientations. -At the lower heat flux levels the 
highest AT Q occurs with vertical surfaces, while at the higher heat 
flux levels, AT Q is highest with the horlzontal-up orientations. AT Q 
is lowest for horizontal -down orientations for all heat flux ranges ex- 
cept at one data point. This discrepancy may be related to the surface 
temperature oscillations which occurred during that run, as may be ob- 
served in Fig. 27. 

These differences observed between ver. cal and horizontal -up 
orientations may be due to variations In the non-boiling convective 
heat transfer arising with these orientations. Convective heat transfer 
is hl^ier with horlzontal-up than vertical orientation. Therefpre, for 
a given heat generation rate a larger heat transfer rate will occur by 
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conduction and convection with the horizontal-up orientation. This 
difference is more significant at the lower heat generation rates and 
results in a lower AT for the horizontal -up orientation. On the other 
hand, the differences in convective heat transferwith orientation do 
not effect the outcome at the higher heat generation rates because as 
heat generation rates are increased the transient periods get smaller, 
and eventually a point is reached where boiling starts even before con- 
vective currents are Initiated. 

As indicated in Figs. 23 and 27, for F113 and respectively, 
the transient surface tenperatures at hic^er heat flux levels behave 
differently with the horizontal down orientation than with the horizontal 
up and vertical up orientations. With the inception of boiling the 
surface temperature drops sharply. However, since the vapor generated 
Is held against the heater surface by buoyant forces the surface tem- 
perature Increases untl 1 the trapped vapor escapes from the edges of 
the heating surface, at which point the surface temperature reduces sharp- 
ly as boiling again takes place When liquid comes into contact with 
the heater surface. The process then repeats Itself. 

This periodic coverage of the surface with vapor slugs Is observed 
in the photographic films, and is in synchronization with the thermal 
recordings. At highest heat flux levels ired here the surface tempera- 
ture drops somewhat with the inception of boiling, but the vapor genera- 
tion rate Is so rapid that the surface becomes continuously covered with 
vapor slugs, causing the surface temperature to Increase sharply toward 
the burnout point, as is noted in fig. 27. 
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The effect of reduced gravity on the surface temperature transients 
is shown in Figs. 21-23 and 25-27 for all three orientations in FI 1 3 and 
LN2» respectively. The level of AT Q is considerably reduced with a 
reduction in gravity for the horizontal up and vertical orientations, 
but is virtually unchanged with the horizontal -down orientation. With 
the inception of boiling at zero gravit} the surface tenperature begins 
to decrease but then the generated vapor bi&bles remain on the surface, 
preventing any further drop in surface teaperature and eventually causing 
the surface temperature to Increase, as shown for example in Fig. 21. 
Since the vapor bubbles are held against the surface by body force with 
the horizontal down orientation at a/g^ reduction of gravity has little 
effect other than the elimination of the oscillation in the surface 
tenperature created by the periodic escaping of vapor slugs at standard 
gravity. 

As Figs. 21 , 23, 25 and 27 illustrate, a continuous increase in 
surface temperature occurs for the horizontal down orientation at both 
a/g a l and a/g°0, while this occurs only at a/g«*0 for the other orien- 
tations, indicating the Invariance of the results with orientation at 
zero gravity. The high speed photographs show that bubbles stay on the 
surface at zero gravity, at least for the period of the observation in 
these experiments. This would tend to reduce heat transfer from the 
surface to the liquid which, for an imposed heat flux, leads to an in- 
crease in surface teoperature, as is noted In Figs. 23 and 27. 

Though not as pronounced as the temperature overshoot, a tempera- 
ture undershoot Is observed in some cases in Figs. 20-26. Different 
speeds of recovery to the steady state values pay be noted. Similar 
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observations of temperature undershoot have been reported [6,7, and 11]. 
At the initial inception of boiling, it is believed that the higher sur- 
face and liquid film superheat activates more sites than the steady 
state boiling can sustain, and the surface temperature reduces below 
the steady state temperature for the given heat flux. The niraber of 
active sites then reduces and the surface temperature rises to the steady 
state level. 

The transient surface tenperatures of two different test surfaces, 
PCG 23 and PCG 29, for the vertical orientation with FI 1 3 are plotted In 
Fig. 20. The temperature overshoot of PCG 29 is higher than the over- 
shoot of PCG 23. This is more apparent in Fig. 29 where the results 
are plotted in the form of AT Q vs. (Q/A)y. The only difference In the 
experimental conditions for these two sets of data was the rou^iness of 
the surfaces used. The pyrex substrate of PCG 29 was polished before 
being coated with the thin gold film whereas the substrate of PCG 23 
was coated as it came. The effect of polishing is indicated by the RMS 
profilometer readings, which varied between 1. 6-2.0 microinches for 
PCG 23 and .09-. 1 microinches for PCG 29. Thus, the surface temperature 
overshoot is larger with smoother surfaces than rougher ones. This is 
reasonable, if the profilometer surface r outness measurements provide 
an indication of cavity sizes influencing the nucleation characteris- 
tics. As will be shown below, smaller cavities require hither surface 
temperatures to become active nucleating sites. 

The transient surface temperature ♦or both and FI 13 are com- 

bined in Fig. 28 for the same surface, PCG 29, with the vertical orien- 
tation. For a given heat flux not only the temperature overshoot but 
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the steady state surface siperheat as well are higher with FI 13 than 
with LN 2 . A comparison of the te^erature overshoot between two fluids 
may also be made In Fig. 29. 

The differences of AT sa £ with FI 1 3 and LN^, and also with PCG 23 
and PCG 29 for FI 1 3 can be explained on the basis of the relationship 
between AT sat , the surface characteristics , and the liquid properties. 
Assuming that a vapor nucleus trapped In a cavity is hemispherical In 
shape and has the same radius as the cavity. It Is possible to calcu- 
late the superheat necessary for the nucleus to grow fro© the cavity, 
using the analysis as presented by Hsu [43]. The pressure difference 
between the vapor within the bubble and the surrounding liquid is given 
by 


?. -n 


a<3‘ 


(13) 


The Clauslus-Clapeyron equation relates the vapor pressure to the phy- 
sical properties as 


AP _ Hfa 
AT " T(l /*-&) 


04) 


Combining Eq. (13) and (14), and assuming that W will give the 

mirlmum superheat necessary for a hemispherical vapor bubble to grow. 


ATut = Te - T»at 


ZCTiat 


05) 


Equation (15) can be used to explain qualitatively the effect of surface 
roughness and the different liquids on AT sat and AT 0 . Since the 
cavity radii are not known, and the surface used In the experiments most 
likely do not contain cavities of a laiforo size* Eq. (15) caimot be 
used for a direct check. 
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Equation (15) shows that, for a given fluid and heater surfaces 
having different roughness, the heater surface superheats necessa»y for 
the Initial formation of vapor bubbles should be Inversely proportional 
to the prevailing cavity sizes, as 


^ATmt ) x _ ( )i 

( ATut ) 2 ( Rc )i 

Assigning that the profllometer RMS values of roughness are Indicators 

for R £ for the two surfaces, and substituting the mean measured RMS 


(16) 


values for R, 


(A Tut)* * - 1.6 

(ATur) u * .03 


20 


07) 


The maaerlcal result given by Eq. (17) can only be Interpreted as an 
Indicator, and has no meaning other than indicating that smoother sur- 
faces require higher superheat for the Inception of boiling. This can 
account for the larger overshoot observed for the surface having the 
smaller RMS profllometer readings, PC6 29, In Fig. 20. 

Equation (15) also can be used to compare ^T sat (Tj B -T sat ) *** 
qulred for Inception of boiling with FI 13 and LNg. The ratio of the 
heater surface superheats necessary for the Initial formation of vapor 
with the two fluids is given from Eq. (15) as 


l ! Tut )m g I <T Tmt \ // (j T14 T \ 

(ATur ^ } pfn / [ J LUt 

Siistituting the appropriate properties into Eq. (18) gives 

ftT...),ii (18.07) («7) Atti )(l«» ... 

{ A Tur) i>lt Ui.lUH-li)/ (W.4XW-*) W 


( 18 ) 


(19) 


Inception of the first dlte is shown by the syrtxrt IS in Fig. 20, end 


was determined from viewing the high speed photogr$>h$. These indict* 
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that boiling Is Initiated before the surface temperature reaches Its 
raaximun value. In other words, (Tj g - T s.t> Is less than AT Q 
^ T max' T sat^‘ Fi 9 ures 30 and 31 are plots of T Ig -T $at vs. heat flux 
for F113 and LNg, respectively, and Illustrate the effects of heat flux, 
orientation, gravity and surface roughness on superheat for Incipient 
boiling. These Indicate that T IB" T sat 1s effected only by surface rough- 
ness and the liquid used, which Is consistent with Eq. (19), which pre- 
dicts that the superheat Is five times higher with FIT 3 than with 
for Inception of boiling. The experiments Indicate that the superheat 
with H13 Is 3-4 tines higher than with U^. Considering the approxi- 
mations made for the derivation of Eq. (15) this Is very reasonable 
agreement. 

In an approach similar to this, a value of 9.4 was calculated for 
the ratio of superheat required for the Inception of tolling In and 
LH 2 (liquid hydrogen) [41]. The experimental values of this ratio varied 
between 2.3 and 9.1, depending upon the surface polish. Polished copper 
surfaces gave results closer to the calculated values than polished 
stainless steel or 600 grit steel surfaces, supporting the results ob- 
tained In the present research. 

2. Delay periods 

In this section, the experimental results related to delay periods 

Tt Tib* ( 7i * 7 V» T3* ( 7i ♦ Tz + T 3 )* percent toiling area 
coverage and toiling spread velocities are presented and discussed In 
chronological order. 
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1. Delay period, fy. 

This period Is the conduction dominated regime. The predicted tem- 
perature rise of the gold film during this period Is plotted In Fig. 32 
and 33 from Eqs. (48) and (49) of Appendix A, for F113 and LN 9 , respec- 
tively, along with the experimental values obtained. The derivation of 
these equations, which are solutions of the problem of two senrl-lnflnlte 
solids with a step Increase In a plane heat source at the boundary of 
these two solids, Is given In Appendix A. The difference between the 
predicted tenperature rise and the experimental value Increases with 
time, the experimental temperature being below the predicted values. The 
linearity of the experimental data and the good comparison at small 
times Indicate that heat transfer Is In the conduction regime during the 
early time periods. 

The next process which occurs In a gravitational field prior to the 
Inception of boiling Is the onset of natural convection. The detailed 
analysis of the theoretical relations predicting the delay time for the 
onset of natural convection as a function of heat flux for the hori- 
zontal up and vertical surfaces are given In Appendix A. 

Equations (54) and (55) In Appendix A are the theoretical relations 
predicting the onset of natural convection for the horizontal up orien- 
tations In and F113, respectively, for a given critical Rayleigh 
meter, Ra»800. The theoretical relations with critical Ra makers of 
1100 and 1700 are tabulated in Table 4 In Appendix A. As described In 
Appendix A, a critical Ra»80G corresponds to a constant wall heat flux 
while Ra*1100 and 1700 correspond to constant wall temperature* with 
free and solid upper boundaries, respectively. Assuming that tin end 
of the conduction, period narks the beginning of natural convection. It Is 



-63- 


possible to compare the experimental results showing Its variation with 
heat flux with those predicted by the equations tabulated In Table 4 In 
Appendix A. These are shown In Figs. 34 and 35 for F1 13 and LN 2 . respec- 
tively, for horizontal up orientations. The agreement appears to be quite 
good. The experimental results in both Figs. 34 and 35 fall between tiie 
two theoretical relations obtained using the critical Rayleigh numbers 
of Ra*800 and Ra a 1100. This Is reasonable, since the boundary condi- 
tions In the tests conducted In the present study are neither constant 
wall temperatures nor constant heat flux In a strict sense. This result 
provides further Indication for the existence of a conduction dominated 
regime for horizontal up orientations. The existence of a convection 
regime prior to boiling has been observed elsewhere [11, 42]. In a 
photographic study with a horizontal up orientation, colusnar Instabi- 
lities rather than the classical Benard cells were observed [11]. Using 
five different liquids, It was established that free convection begins 
after the critical Rayleigh nusber of 1700 Is surpassed [42], 

Equations (57) and (58) In Appendix A are the theoretical relations 
predicting the onset of natural convection for vertical orientations In 
LN 2 and F113, respectively. Figures 36 and 37 show the comparisons 
between the experimental results obtained with vertical orientations In 
F113 and U 2 with those predicted by the theoretical equations ( 57 ) 
and (58), respectively. In both cases, the agreement Is reasonably good 
for the lower levels of heat flux. As heat flux is increased, the experi- 
mental values deviated from the theoretical predictions more and more. 

For both liquids , the experimental values were always lower titan the 
predicted results, and may be anticipated on the basis of the following; 

As explained In Appeal x A, the transition from conduction to convection 
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on a vertical surface begins only when some effect from the leading edge 
Is propagated up the plate to the particular point in question. Before 
this leading edge effect has propagated to that point, the region ef- 
fectively does not know that the plate has a leading edge. Therefore, 
the delay time f-\ In Eqs. (58) and (59) Is the total time necessary for 
the convection to be established over the entire surface. It Is the 
elapsed time necessary for the disturbance coming from the leading edge 
to reach the upper edge of the surface. Convection at a point begins 
when the disturbance coming from the leading edge reaches thlt point, 
and with the transient heating process taking place, a reduction In the 
rate of Increase of the surface temperature nay take olace before the 
delay time has reached the maximum possible value, corresponding to the 
vertical height of the surface, due U the Increased heat flux being 
carried away by convection. The experimental values pressed In Figs. 
36 and 37 correspond to that time when the rate of Increase In the 
average surface temperature begin to reduce, aid the experimental values 
of T] are less than predicted, as expected. In addition, as the heat 
flux level is increased, boiling nay start before convection has become 
established over the entire surface. When heat fluxes are sufficiently 
high, the surface temperature may rise so rapidly that boiling Is es- 
tablished before any convection can take place. This Is probably why 
the existence of a convection dominated regime was not observed In some 
experimental works [2,8]. Heat flux levels In the range of 10 5 -2xl0 6 
3TU/HR-FT 2 were used In experl amts with water [2,8], which goes above 
the u urn out heat flux, whereas the heat flux levels were kept below the 
burnout point In the present experiments and varied between 3x10^-1.5x10 
BTU/HR-FT 2 with and P113. 


Jh> 
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11, Delay period for Inception of boiling, 7 ^: 

T^g Is the time elapsed from a step change In heat flux to the 
appearance of the first bubble. Figures 38 and 39 Illustrate tht effect 
of heat flux and gravity on the elapsed period, 7 j a , for vertical sur- 
faces In FI 1 3 and LN^ , respectively. The results indicate that 7 ^ 
decreases with increasing heat flux and with reduction of gravity. The 
difference In the 7 jg values at a/g»l and a/g °0 become less as heat 
flux Increases. One might expect that at sufficiently high heat fluxes 
the body forces no longer will have an effect. These effects are re- 
lated to the amount of energy available to heat the liquid film adja- 
cent to the test surface and can be explained as follows: 

Hsu and Graham [42, 43] have postulated and experimentally veri- 
fied that a bubble will t,»: grow until the liquid surrounding the bubble 
nucleus J .i wanner than the equilibrium superheat temperature as given 
by Eq. (15). Thus, the delay time ^IB depends on the transient growth 
of a • '.;y thin thermal layer adjacent to the heat 4 ng surface. As heat 
flux Increases, the liquid film tenperatur-- Increases mure rapidly and 
T ib becomes smaller. Similarly, with the reduction of gravity con- 
vection effects are reduced or eliminated, and more energy Is available 
to heat the liquid film adjacent to the surface. As heat flux Increases, 
the delay time shortens and the effect of convection is reduced, and the 
effect of gravity on Tj B this becomes smaller. 

The values of ^IB for a vertical surface In FI 1 3 and at a/g*! 

and a/g*0 are presented In Fig. 40 along with the de*ay times "or the 
onset of spreading of nucleate boiling. The delay time, ^IB’ Is con- 
siderably shorter with than with Fill. Although the therwl differ 
slvltles of both liquids are similar* FI 13 stores * heater surface 
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superheat about five times larger t an LN ? for the inception of boiling, 
accounting for the longer delay time with F113. 

The flagged symbols in Fig. 40 are the delay times from the step 
inDut in power to the beginning of spreading of the active boiling sites, 
as observed in the high speed photoaraphs. Boiling initiates first at a 
few sites, and the increase in the n'mber of active sites on the surface 
does not take place until an additional delay period. These additional 
delay periods are the differences between the flagged and unflagged sym- 
bols in Fig. 40. The surface superheat, when spreading begins, is 
higher than the incipience of the first sites. This phenomena occured 
only with LN^. With F113, the spread of boiling sites immediately 
followed the appearance of the first bubble site. It is also 'served 
in the photogriphs that the inception of boiling on vertical surfaces 
begin at t? . per edge and penetrates downward. This can be seen in the 
film sanples shown in Fig. 95 of Appendix C. As discussed earlier, con- 
vection occurs at the upper part of the test surface at a later time, 
which allows the liquid boundary layer near the upper edge to be heated 
to a higher temperature than that below. In addition the lie - i film 
tenperature in the upper pare Is always higher due to the upward motion 
of the convection currents. Therefore, boiling may be ejected to begin 
preferentially at the upper side of the surface provided that the cavity 
size distribution is uniform over the entire surface. On horizontal 
surfaces, boiling is observed to becr'n near the center of the test sur- 
faces and spread outward. Vest and Sweeney [44] showed, icing holography, 
that the liquid tei^erature Is hipest at central lccatlons for a heated 
horizontal surface facing upward, which explains the fact that the Ini- 
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tiation of first boiling sites takes place at the central part of the 
test surfaces facing upward. 

iii. Delay period for maximum surface temperature, { Ti ♦ T 2 ) 

The test surface temperature rises at Its highest rate during the 
conduction period, and tapers off with the onset of convection. It 
tapers off further with the activation of the first few nucleating sites, 
but still continues to increase until a certain percent of the area is 
covered with active boiling sites. 7jg is the delay time from the 
power input to the activation of first site, and ( T-| + T 2 ) is defined 
as the delay time from the power Input to the poi..c where surface tem- 
perature reaches its maximum value. These are compared for various levels 
of heat flux on Figs. 41 and 42 for a vertical surface with FIT 3 at 
standard and zero gravity, respectively. These indicate that the surface 
temperature continues to rise even after the activation of the first 
nucleating site. This is also shown in Fig. 20 where activation of the 
first site is identified by the syfd>ol “IB". Actually, as shown in 
Figs. 95 and 96 in Appendix C, the test surface temperature continues to 
Increase unti 1 enough sites are activated such that the heat transfer 
by conduction, convection and boiling is equal to the power input. 

( T| + T 2 ) depends on: (1) the rate of Increase of temperature on 

the nonboiling portions of the heat transfer surface and (2) the rate 
of spread* ng of the active sites. The surface temperature as measured 
is a mean value, and It is determined from the total .-esi stance of the 
surface. The surface resistance at any time is equal to the sun of the 
resistance of the boiling part ftg and u*e resistance of the nonboiling 
part R^g of the surface. A relationship between the boiling and non- 
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boiling surface temperatures and the boiling and nonboiling surfac 
areas can be found for the condition when the surface temperature reaches 
its maximum value. The relation is given by Eq. (81) in Appendix E, 
which indicates that the surface temperature continues to Increase until 
a certain percent of the area is covered by active sites. 

Figures 43 and 44 show the effects of heat flux and gravity on 
( 7*-] + Tg) for surfaces in FI 1 3 In the horizontal up and horizontal down 
oriental jns , respectively. As Figs. 41-44 illustrate, ( 7] + Tg) de- 
creases a (Q/A)y increases. The difference between the energy gen- 
erated and the heat transfer by conduction and ..jnboiling convection 
becomes larger as (Q/A)y is increased. This In turn increases surface 
temperature at a faster rate and b.ings the liquid film superheat to a 
boiling level sooner. 

Figure 43 and comparisons of Figs. 41 and 42 show that reduction in 
gravity reduces ( T-| + Tg) for the vertical and horizontal up surfaces. 
Comparing Figs. 41-44 also shows that ( 7] + Tg) shortest for the 
horizontal down orientation. The reduction of gravity has an effect 
similar to that of rotating the surface to the horizontal down orien- 
tation in that the convection part of the heat transfer is reduced or 
eliminated. Therefore, as may be seen in Fig. 44, ( Tj ♦ Tg) does not 
vary slgilficantly from a/g=l to a/g°0 with the horizontal down orien- 
tation. The explanation for the effect of gravity and orientation on 
( 7”i + T 2 ) is similar to that given for *7^, and is due to the influ- 
ence in convection effects. 

Behaviors similar to these occurred with although the values 
of ( 7} + 7 2 ) are considerably shorter. These are shown in Figs. 51-53, 
along *1th other delay times to be discussed shortly. 
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iv. Spreading periods, T 3 and X s ‘ 

'Tj is the elapsed time from the point when surface temperature 
reaches its maximum value to the time when spreading of nucleate boiling 
over the surface ends. T s is the elapsed time from the beginning to the 
end of spreading of nucleate boiling over the surface. For most cases 
the difference is negligibly small. Since, as explained before, surface 
temperature continues to Increase after spreading started is always 
greater than ~Cy 

The starting point, direction and velocity of boiling spread vary 
depending on heat flux, orientation, liquid media and the test surfaces 
used. Figures 45 and 46 are plots of heat flux versus elapsed tine, 'X y 
for F113 with vertical and horizontal up orientations, respectively. T 3 
values for LN 2 for various orientations are shown in Figs. 51-53. They 
show the effect of gravity and (Q/A) T on'Vy which decreases with Increas- 
ing (Q/A) T and reduction of gravity. Increasing (Q/A)y Increases the 
rate Ci heating of the boundary layer, producing more superheat which, 
in turn, speeds the activation of the sites on the surface. Reduction 
of a/g decreases the convection heat transfer, which also Increases the 
rate of heating the boundary layer of sites, and thus also reduces T 3 . 

Figures 59-61 and 62-64 show the effect of heat flux and gravity on 
T s for various orientations in F113 and LMg. respectively, in terras of 
the relationship between time and the fraction of surface area covered oy 
nucleate boiling, as observed from the high speed photographs. These will 
be considered further shortly. The effects of heat flux, gravity and 
orientation on T s in LH 2 are similar to the effects in F113. The differ- 
ence Is that values of T s are smaller with U 2 than F113, i.e. fear ver- 
tical orientation with (Q/A) T »2800 BTU/HA. : Z^)ni3/(7^)Ui2*.13/.03M.33. 
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As Eq. (18) Indicates for a gl ven cavity s1ze» about one fifth the 
superheat Is needed In LN 2 as compared to F113 for a bubble to grow. This 
may be the reason for the faster spread of active sites In than In F113. 

v. Delay period ( 7 ^ + 7 ^ + 7 * 3 ) : 

For some thermal systems, particularly nuclear reactors, the elapsed 
time between the Initiation of power input and the complete coverage of 
the test surface by boiling sites could be very critical. This period 
Is represented by 7 ^ ♦ 7 ^ + 7 ^. Figures 47-49 show the effect of 
heat flux, gravity and surface toughness on the total elapsed time, 

7*1 + 7 ^ + r 3 . f° r Mch orientation with F113, while Fig. 50 Is a 
composite plot of the results pres hi ted In Figs. 47-49. Figures 47 and 
48 Indicate that the total delay time decreases with Increasing heat flux 
and with reduction of gravity for the vertical and horizontal up orien- 
tations. While reduction of gravity has a considerable effect on the 
total delay time with the horizontal up and vertical orient tlons. Its 
effect with the horizontal down orientation Is Insignificant In Fig. 49. 

This phenomena Is more apparent In Fig. 50, which also Indicates that at 
a/g»0 the total delay time Is Insensitive to orientation. Figures 47- 
50 also Indicate that a reduction of gravity has less and less Influence 
as heat flux Is Increased. Since the total delay time 7 i + T2 + T3 
Is the stsn of the Individual cocponents 7 i» Tg T3* 
tlon for the effects of heat flux, gravity and orientation on the total 
delay tire Is the same of that given for the Individual components, and 
Is due primarily to the Increased rate of heating of the boundary layer 
resulting either from the Increased (Q/A) T or the decreased nontolling 
convection effects. 
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Experlmental results of the total delay tire obtained with LN^ are 
presented in Figs. 51-53 for each orientation at a/g a l ami a/g^O. These 
include plots of ( 7"i + 7^ an ^ T 3 ^ or ^2’ were d ^cussed earlier 
in connection with FI 13. The general results are similar to effects ob- 
served with FI 1 3 . 

The total delay ttme results obtained with LN^ and FT T 3 at a/g=l 
for each orientation are presented in Figs. 54-56 for comparison. Re- 
sults obtained at a/g=0 for all orientations are presented together in 
Fig. 57. These figures indicate that the total delay time in FI 1 3 is 
about 5-6 times that for L 1 ^ at any given heat flux. The difference in 
the total delay time in and FI 13 is most likely due to the property 
diff rences of surface tension, latent heat and saturation temperature 
which gave rise to inferences In ( 7^ + 7^) and 7^. 

Both ( 7^ + T^) and 7*3 are larger for FI 13 than for Lf^. 

A rapid decrease in the total delay time with increasing heat gen- 
eration rate was also reported for water [8] , and is plotted with the 
results of the present study In Fig. 58. Comparisons between these are 
tenuous due ♦o the differences In the test surfaces and liquids, but the 
linearity of log-log scale indicating exponential relations is Inter- 
esting, and perhaps can be exploited for estimating purposes. 

3. Percent area coverage and spread velocity 

Figures 59-61 show the spreading of active sites over the test sur- 
face as the percent area covered versus time, at a/g°l and a/g=0 for three 
orientations in FU3. As heat flux is Increased and gravity Is reduced 
the boiling sites activate at a faster rate and the test surfaces are 
covered by active sites In shorter time. Comparison of Figs. 59-61 in- 
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dicate that a reduction of gravity has the highest effect In reducing 
the area coverage time with hor^pntal up orientation and lowest effect 
with horizontal down orientation, with a reduction In this effect occurlng 
as heat flux increases. This Is consistent with the relationship between 
natural convection and 7*3 discussec previously. 

The percent area coverage by active sites along with the delay time 
for Inception of boiling in at a/g»l and a/g=0 for the vertical and 
horizontal up orientations are shown in Figs. 62 and 63, respectively. 
Percent area coverage for the horizontal down orientation Is shown in 
Fig. 64. Increasing the heat flux has similar effects on the spreading 
of active sites In LN^ and FI 1 3 , In that T 5 decreases. The effect of 
gravity on the boiling spread Is more pronounced with FI 1 3 than with LN,,. 
In fact, upon close examination of Figs. 62-64, one can say that reduc- 
tion of gravity has no effect on either the Incipient boiling point or 
the rate of formation of active nucleating sites. If this Is a general 
characteristic of cryogenic 11qui<fe, this could have significance In 
future research directions connected with heat transfer for space appli- 
cations. After the inception of first site boiling spreads 4 : 'Jl 2 at a 
faster rate than in F113. The effect of orientation on boiling spread 
is more pronounced with F1 13 than with The explanation for the 

effect of (Q/A) gravity and orientation on the boiling spread is similar 
to that given for TyWd is due to increased heating rate of the boun- 
dary 1 tyer by either the increased heat flux or reduced convection ef- 
fects. 

Comparisort in the rate of spreading with the two test surfaces PCG 
23 and 29 In FU3 may be made In Fig. 59 for the vertical orientc ion. 
Conparison between PCG 24 and 29 with the horizontal down orientation 
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may be made in Fig. 61. Profilometer RMS values for PCG 29, 23, and 24 
are between .09-. 1, 1. 6-2.0 and 3.0-10.0 raicroinches, respectively. 

As shown In both Figs. 59 and 61, after .he inception of the first site, 
additional boiling sites activate or spread at a faster rate on the 
smoother suraces. This is as might be expected. According to Eq. (15), 
smoother surfaces, which are presumed to have smaller cavities, require 
higher superheats for the inception of boiling, thus raising the test 
surface temperature to a higher value, which in turn speeds the spread 
once boiling is initiated. 

The effect of heat flux and gravity on 7 ^, ( , 

T'i + T2 + T3 and percent area for a vertical surface in FI 1 3 are 
shown together for convenience in Fig. 65. The arrows on Fig. 65 indi- 
cate the percent of the surface area covered by the active sites at the 
delay periods fry and ( 7 ^ + . The percent area covered by active 

sites at the delay time ( 7 ^ + 7 ^) increases with heat flux and with 
reduction of gravity. Up to elapsed time ( f-[ + 7 * 2 ) surface tem- 
perature Increases and after this time the surface temperature decreases. 
The slope of the curves is a measure of the rate of boiling spread, 
which Increases with increasing heat flux and with a reduction in gravity. 

The average velocity at wnich the spreading of nucleate boiling 
took place was determined from the high speed photographs. Two direc- 
tions are chosen from the point at which the Initial bubble forms, one 
cor^. ponding to the maximum rite of spreading, the other the rainitam 
rate of spreading. The spreading speed in these two directions are then 
averaged. The spreading speed In any direction is not constant but 
varies from region to region , giving a large deviation frvn average 
values. Figure 66 shows the average boiling spread velocity for all 
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orlentatlons In F113, at a/g** 1 and a/g=0, and Fig. 67 shews the corres- 
ponding results for LN^. These can be used for qualitative purpose only, 
however, because of the relatively large variations observed, whl Ji &r° 
represented, typically, by the vertical dotted lines In Fig. 66. The 
spread velocity Is faster with LN^ „han FI 1 3 and, for a given heat flux, 
the spread velocity Increases with reduction of gravity. In FI 9s. 66 and 
67, the darkened symbols correspond to a/g=0 for the same orientation, 
at a/g=*l represented by the open symbols. In all cases, both fluids and 
all orientations, the darkened data points lie above the open ones, in- 
dicating that the spread velocity Increases with reduced body force. 

The explanation for this effect was given before with the discussions of 
T 3 .7 ' s and the percent area coverage. Since the spread velocity is re- 
lated to > which is the boiling spread time, the same factors govern 
both. 

B. STEADY STATE NUCLEATE BOILING HEAT TRANSFER RESULTS 

In the steady state part of the present study the following para- 
meters were measured: 1) heat flux, 2) test surface temperature, 

3) active site and population density, 4) bii>ble frt ,uency, and S) 
maximum bubble si2e. The step power Input, body force (gravity) and the 
liquid used were independent variables. The results are discussed under 
the related headings. Comparisons between some of the existing nucleate 
boiling models with results of the present work are also Included In 
this section. 
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1. Heat Transfer Results 

The heat transfer results are presented In the form of plots of 
heat flux versus superheat for different orientations and body forces. 
Figures 68-70 show the heat transfer results obtained at a/g=*l and a/g a 0 
for the three orientations in F 1 13. They also include the natural con- 
vection correlation equations for comparison purposes. The effect of 
the variou. orientations are compared in Fig. 71 for a/g=l and in Fig. 72 
for a/g=0 Figure 73 is heat flux versus AT $at for LN ? and illus- 
trates the effects of both orientation and body force on nucleate boil- 
ing heat transfer. 

Experimental results obtained with two different surfaces, PC6 23 
and 29, at about two months apart are presented together in Fig. 68 for 
the vertical orientation and in Fig. 70 for horizontal down orientations. 
The consistency of these results, which were obtained with different 
surfaces at different times, demonstrate the reproducibility of the 
data. 

Figure 71 Indicates that, for a given AT at lower levels of heat 
flux, the horizontal down orientations provide the highest heat flux, 
the horizontal up provide the lowest, and the vertical orientation falls 
between, for F113 at a/g»l. Similar results were obtained with and 
are presented in Fig. 73. These results showing the effect of orienta- 
tion are consistent with the results of Refs. £ 24 - 28 ], which also show 
that results obtained with the thin gold films compare well with those 
obtained with conventional metal heat transfer surfaces such as ribbons 
and discs 7 [ 24 - 28 ]. Figure 2 of Ref. [ 25 ], in which a copper disc was 
used as test surface. Is replotted in Fig. 74 for comparison. The var 
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latlon in heat flux with orientation for a given AT is believed re- 
lated to the effect of “microlayer evaporation* at the base of the babies. 
At the lower levels of heat flux, with the horizontal down orientation 
buoyant forces tend to hold the vapor bubbles closer to the surface, thus 
increasing the size of the nricrolayt*. With the horizontal up orien- 
tation the buoyant forces tend to decrease the size of the microlayer 
trapped beneath the vapor bubble, thus decreasing the heat flux for a 
given heater surface superheat. This is consistent wi th the experimental 
observation by Hospetl and Nesler [32] that heat transfer associated 
with a hemispherical bubble Is greater than that with a spherical bubble, 
due to the larger microlayer trapped with hemispherical biibles than 
the spherical ones. Di>° to the buoyant forces, the horizontal down or- 
ientation tends to produce hemispherical bubbles while the horizontal 
up position favors spherical bubbles. 

Experimental results for all three orientations In F113 and LN 2 
consistently fall above natural convection correlations, Figs. 68-71 , 73. 
Photographic results Indicate that boiling was present for all data 
points shown in these figure, and since the heat flux is higher with 
nucleate boiling than with natural convection, this Is to be expected. 

The lowest data points in Figs. 68 and 70 are results obtained with in- 
duced boiling as indicated. The Induced boiling points correspond to 
such low heat flux levels that a step increase to that level does not 
result In nucleate boiling on the test surface, but once boiling has 
been Initiated at a higher heat flux level, It continues when the heat 
flux Is reduced to that level. This procedure for obtaining boiling at 
lower heat flex levels Is an example uf the hysteresis effect, which 
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moves the boiling curves toward left (e.g. Ref. [41j). Figure 71 indi- 
cates that the horizontal down orientation requires the lowest heat flux 
to sustain the induced boiling and the horizontal up requires the high- 
est. This is related to the effect of natural convection heat transfer, 
discussed In Section A-2 for the transient process. 

As Indicated In Figs. 68-70, for a given heat flux the boiling sur- 
face superheat Increases with the reduction of gravity from a/g= 1 to a/g°0 
for the vertical and horizontal down orientations, whereas it decreases 
with the horizontal up orientation. Reduction of gravity reduces or 
eliminates the natural convection effects, and It might be expected that 
the nucleate boiling performance should be Independent of orientation. 

That this Is the case Is Illustrated In Fig. 72 for F113, where heat flux 
versus surface superheat is plotted for all orientations at a/g°0. 

2. Test Surface Superheat Results 

Figures 75-79 are plots of surface superheat versus time for the 
different constant levels of power Input and the various orientations 
used with FI 1 3 and LN^, and show the transient surface temperatures 
during the transition from steaqy state at a/g=l to the lower gravity 
level, a/g»0. It is possible that a new true steady state consition Is 
not yet reached during the 1.4 second free fall period, but examination 
of the surface temperature behavior in any of the Figs. 75-77, for 
exauple, would appear to support the idea that a new steady state Is 
reached even within such a short tick. The experimental results of sur- 
face temperatures at a/g*0 presented In Figs. 68-73 were obtained from such 
plots as Figs. 75- 79, and show that surface superheat increases with a 
reduction of gravity for vertical and horizontal down orientations (see Figs. 
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75, 76). The effect of reducing the gravity drastically becomes more 
pronounced as the heat flux Is Increased. As Illustrated In Fig. 78, 
when the heat flux Is high enough at a/g*0, the surface temperature In- 
creases continuously toward the burnout point. Similar results were re- 
ported by Littles [24] in experiment- with F113 at atmospheric pressure, 

In that an Increase in surface temperature occured with the vertical and 
horizontal down orientations, and a decrease occured with the horizontal 
up orientations as gravity was reduced from a/g"l to a/g**0.01. Figure 77, 
where surface superheat versus time for the horizontal down orientation 
Is given. Illustrates that the surface tenperature fluctuations caused 
by escaping vapor bubbles from the sides at a/g»l are danped by a reduc- 
tion of gravity to a/g^O, and reappear when standard gravity condition 
Is resumed after the impact. The different behavior of the test surface 
temperatures with reduction of gravity for the various orientations are 
caused by combinations of several phenomena: 1) reduction of gravity 

reduces or eliminates the convection heat transfer, Increasing the ’•ate 
of energy transfer by boiling alone which, In turn, increases r • AT $at 
necessary. 2) As shown In Fig. 98 in Appendix C, which Is a typical 
saople of high speed photographs of steady state boiling at a/g=0, 
departure of vapor bubbles from the vicinity of the heating surfaces ceases 
shortly. after gravity Is reduced from a/g=l to a/g=0, and the bubbles 
remain on the surface and continue to grow there. Vapor bubbles remain- 
ing In the vicinity of the heat transfer surface behave as thermal Insu- 
lators, and this the test surface tenperature Increases. 3) The third 
phenomena may be a succeeding consequence of the first two effects. As 
the surface temperature increases, sites 'with smaller critical radii 
become activated. The addition of new active sites then decreases the 
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average test surface temperature, as was discussed earlier wits regard 
to transient boiling. 4) The reduction of gravity from a/g*l to a/g a 0 
makes the b<i>ble shapes similar for the different orientations, during 
the growth period, because of the elimination of the buoyant forces. 

The bubble shapes become .»re hemispherical for all orientations at 
a/g“0. Thus, the ir cro’ayer trapped under the growing vapor bubbles 
Increases with a reduction of gravity where this hemisph *rical shape Is 
promoted. The maximum Increase occurs with the horizontal up orienta- 
tion, where bubbles growing at a/g°l are primarily spherical. The In- 
creased mfci-olaye*' area means more rapid llqu) ; s>,aporat1on, and a hlghe- 
rjte of heat transfer resulting in a largt .r^ace tenperature drop. 

The net results of these opposing phenomena depends on the initial orien- 
tation at a/g«*l. 

3. Active Site Density and Bubble population Results 

Figures 30 and 81 present the results of the pnotographlc study In 
which the active site densities in nucleate boiling were determined for 
FI 1 3 and respectively. The results are presented in the form of 
plots of active siu. (tensity versus heat flux for a 79 ® 1 and a/g»0 and 
for the vertical, horizontal up, and horizontal down orientations. The 
active site densities show., represent the number of nucleatlon sites 
pe* unit area which were observed to be active during the entire length 
of the film. Inducted for reference on F-ig. 80 are the approximate 
burnout, minimum incipient and Induced boiling Heat flux values. Figures 
81 and 82 present the average population densities along with the active 
site (tensity for IN,, and F113, respectively. Tne average. populc.Ho>* den- 
sities represent the average of the bidtele popu'at'ons on the suffice at 


any given instant of time. The vertical bars on each data point indi- 
cate the observed variations associated with each of the mean values. 

It is noted t-iat the active site densities are always larger than the 
average population density. This occurs because of the fact that at 
any instant of time a nunfcer of the active sites may be in the "waiting" 
period. 

The results of present investigation shown In Figs. 80-82 indicate 
that the active and population site densities increase considerably with 
heat flux. The linearity of the plots indicates the existence of an 
Aponential relation between the active site density and heat flux such 
as 

<3 

-2sL = Cow*T4MT f ~J — ^ tS» i (20) 

The relationship between heat flux and active nucleation site density has 
been noted by numerous investigators, observing that increasing the heat 
flux results in the activation of additional nucleating sites in order 
to facilitate the transfer of the additional heat. The results of West- 
water and Kirby [45] for carbon tetrachloride boiling on an artificially 
roughened oxide coated glass surface and the results of Judd and Merte [34] 
for F113 boiling on a smooth oxide coated glass surface are shown In Fig. 
83, along with the results of the present work, for purposes of conparison. 
It is noted that a considerable dlscrepancy^exists between the present 
work and that of Judd and Merte [34] (for a given heat flux Ref. [34] 
reports an active site density four times as large) even though a saooth 
glass surface was used In both cases. This discrepancy could reflect 
the difference In the coatings used, a vacuum deposited gold film In the 
present case as compared to a fired oxide coating in [34]. The discrep- 
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ancy could also be due to procedural differences used in determining the 
site (tensities. Figure 97 is a representative series of high speed photos 
taken during the course of the present research for steady state boiling 
on a vertical surface. These show that the distribution of sites Is not 
uniform over the surface. Therefore, depending upon the section of the 
area chosen for site counting, the results could vary manyfold. In the 
present work, as explained before, measurements were made over 50 - 802 of 
the surface, whereas in Ref. [34] wily approximately 32 of the test sur- 
face was viewed for determination of the site densities, it is possible 
that the sample size used in Ref. [34] was not large enough to provide 
a meaningful result for the site density. 

In addition to the effects described in the above, Figs. 81-82 
indicate that even though the change is not large, the active site den- 
sity increases slightly with a reduction of gravity, while the population 
density decreases slightly with a reduction of gravity. The increase of 
N/A is to be expected, if the heat transfer per site remains constant, 
or if it does not increase enou#i to compensate for the increase of heat 
transfer by boiling, owing to the elimination of convection heat trans- 
fer accompanying the reduction of gravity from a/g=l to a/g=0. On the 
other hand, the decrease in population density could result from the 
Inaccuracy in the determination of n/A at 3 / 9 = 0 . Since the bubbles stay 
and grow on the surface at a/g»0, neighboring sites may become covered 
by the adjacent growing bubbles, preventing an accurate count of the 
sites active at any Instant of time. 

Figures 80 and 81 Indicate that orientation, like gravity, do not 
appear to have large effects on the active site densities. The active 
site density is sli^tly larger with the horizontal up orientation and 
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slightly lower with the vertical orientation. This could be related to 
the microlayer effect. Bubbles generated on surfaces at the horizontal 
up orientation generally are more apt to be spherical in shape, and thus 
have a minimum microlayer contribution, and each bubble results in less 
heat transfer. Therefore, for a given heat flux more active sites may be 
required with the horizontal up orientation. On the other hand, as noted 
in Fig. 80, the active site density for the horizontal down orientation 
is slightly higher than for the vertical orientation, which seems contra- 
dictory to the argument given above. Since bubbles generated on the 
horizontal down surfaces are more apt to be hemispherical in shape, and 
thus have larger microlayer contribution, each bubble should contribute 
more to the heat transfer. For a given heat flux fewer active sites 
should therefore be required with the horizontal down orientation. On 
the other hand, for a given heat flux a reduction of convective heat 
transfer increases the rate of energy transfer by boiling. It is 
therefore possible to increase the active site density, provided that the 
heat transfer per site and the frequency of departure remain constant or 
not increase accordingly. Finally, growing vapor bubbles remain close 
to the surface in the horizontal down orientation until they escape from 
sides, and behave as thermal insulators. The surface temperature thus 
increases which, in turn, activates additional sites at the smaller 
cavi ties. 

4. Maxima Bt&ble Size Results 

Bubbles attain their maximum sizes at departure from the surface 
in nucleate boiling with a saturated liquid, whereas in subcooled boiling 
the departure sizes may be smaller than the maximum size attained during 
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the growth period doe to partial collapse before departure from the 
surface. 

Figure 84 is a plot of the maximum bubble diameter as a function of 
heat flux for the horizontal up and vertical orientations with FI 1 3 at 
atmospheric pressures. The vertical bars represent the observed varia- 
tions of the maximum bubble sizes. The results Indicate a slight de- 
crease in average maximum bubble diameters with increasing heat flux. 
Beckman, Merte [46] and Gaertner [47] have also observed a small decrease 
in bubble diameters with increasing heat flux in their experiments with 
water. The decrease in oubble diameter could be due to momentum effects 
associated with the bubble growth. 

Figure 84 also indicates that the average maximum bi&ble diameters 
are larger with the horizontal up orientations than with vertical, in 
FI 1 3, and may be a consequence of the Induced velocity sweeping the bubbles 
off sower with the vertical surface. 

Figure 85 presents results obtained for the vertical and horizontal 
down orientations with LN^. With the horizontal down orientation there 
is no departure diameter In the conventional sense, because buoyant forces 
hold the bubbles against the surface imtil the vapor slugs that cover 
the surface escape from the sides. In the present experiments with LN^, 
the test surface is tilted about 3-5° from the horizontal to help vapor 
bubbles to move off the site more frequently, which then made it possible 
to observe the maximum bubble diameters presented In Fig. 85. The max- 
imum bubble diameters decrease slightly with heat flux In LJ^ as In FI 13. 

As shown in Fig. 85 the average maximum bubble sizes are higher with the 
horizontal down orientations than with vertical. This could be due to 
the distortion of the bubble shapes by buoyant forces, or also for the 
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reason that the bubbles may be swept away sooner by the convection 
currents on vertical surfaces. 

The bubble departure diameter was first formulated by Fritz [48], 
based on the work of Bashfort and Adams [49] from a static force balance 
of surface tension and buoyancy for.®s, who obtained the following eq- 
uation 

D = 0.0119 0 (21) 

Equation (21) is plotted in Fig. 84 for an arbitrary constant contact 
angle 0 n 24° for comparison. Equation (21) has been shown to predict ■ 
departure diameters reasonably well at low heat fluxes [46]. Since the 
dynamic forces were ignored, Fritz' equation does not predict the de- 
parture diameter well at higher heat flux levels. Among others, Stan- 
iszewski [50] found that the bii>ble departure diameter depends on the 
growth velocity in the last stages of growth and suggested a correction 
to Include the dynamic effects, resulting In 

v ‘ ooov ^i^)( l ' 0i35 £) (22) 

Both Eqs . (21) and (22) Include the C /( ft -ft) term as a multiplier. 
This term can be considered to represent the ratio of the surface ten- 
sion to buoyancy forces, and can be used to obtain an indication of the 
expected maximum bubble sizes In two different liquids. Taking the ratio 
of Eq. (21) for two liquids and assuming that the contact angle is the 
same for both 

Dlmz (Tlm - phi 

~ < J> " 3 (fl-ft) LNJ 


Dm3 


(23) 
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SUbstltutlng the properties of F113 and Lf^, In Eq. (23) gives 

122 (24) 

Dfiis 

In saturated boiling the bubbles attain their maximum sizes at de- 
parture, and the experimental results could be used for comparison with 
the value of Eq. (24). substituting the departure sizes fro® Figs. 84 
and 85 for tyA«*9000 BTU/Hr-Ft 2 will give 


26 j 

Dm3 ^ 


(25) 


The result given by Eq. (25) agrees reasonably well with the pre- 
dicted results of Eq. (24). 


5. Bubble frequency results 

Figures 86 and 87 present the results of the analysis of the photo- 
graphs In which the average frequency of departure of bubbles In nucleate 
boiling were determined for FI 1 3 and respectively. The results are 
presented In the form of plots of average frequency of bubble departure 
versus heat flux for each orientation at a/g»l. The results 
Indicate that the frequency of vapor bubble emission varies as a power 
of the heat flux. 

As the photographic samples In Appendix C and Fig. 18 Indicate, a 
bubble period consists of two parts: the active period, or elapsed time 
from the appearance to the departure of a bubble, and the waiting period, 
or elapsed tine between the departure to the appearance of the next bubble 
at the same site. The waiting period depends on how rapidly the new 
liquid boundary layer can be heated to the minimum superheat necessary 
for the next bubble to grow. The required superheat for the formation of a 
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bUbble, as discussed earlier, depends on the cavity size, surface ten- 
sion, saturation temperature, etc. Therefore, for a given liquid the 
waiting period depends on the amount of the excess heat flux over the 
amount transfered by conduction and convection. The growth period of a 
babble depends on the bi&ble growth rate, which Is dependent upon the 
surface temperature and the temperature of the fluid adjacent to the test 
surface, both of which Increase with Increasing heat flux. Therefore, 
as heat flux Is Increased both the waiting period and the growth period 
decrease thus increasing the frequency of bubble departure. As ob- 
served previously, however, the maximum size decreases somewhat. 

Figures 86 and 87 also indicate that the average frequency of bii>ble 
departure is filler with the vertical orientations than the horizontal ip 
orientations. As shown In the previous section dealing with departure 
sizes (e.g. Figs. 84 and 85) the vapor bubble sizes are larger on depar- 
ture for the horizontal down and horizontal up orientations than for the 
vertical orientations. If the growth rates remain coeparable for the 
different orientations, then the growth period will be longer for the 
horizontal surfaces than for the vertical surface, which will In turn 
reduce the frequency of depa»*ture for the horizontal surfaces. In addi- 
tion, buoyant forces tend to hold bubbles on the surface with the hori- 
zontal down orientations, Inhibiting their departure, and thus become the 
major contributing factor In reducing the frequency for that orientation. 
Furthermore, for a given heat flux the energy transfer rate by boiling 
Is less with the horizontal up orientation than with the vertical orien- 
tation due to the higher convection heat transfer which. In turn, results 
In a reduction in the frequency of departure. 
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Figure 86 Includes the results of [34 j for F113 for the sake of 
comparison. Considering that the present results were obtained at about 
2° subcooling, which, according to the results of Ref. [34] reduces the 
frequency of bii)ble departure, agreement with [34] Is excellent. 

As reported by Jakob [51], when the frequency of bubble departure 
was plotted against the departure diameter a hyperbolic relation was 
obse:*ved such that 



CoMJTANT 


(26) 


Later Zuber [52] proposed the following relation Involving surface ten- 
sion and buoyancy forces: 


•T. “D — 1 16 ^CT ( fj ~ Pv ). 
T ' 2 


A 


s 


(27) 


McFadden and Grassmann [53] demonstrated the Inadequacy of Eq. (27) for 
their experimental results with U^. They developed the following rela- 
tion with dimensional analysis, coupled with available experimental data 


f * * C'JHSTAM t (28) 

Later Stanlszewskl [50] and Perkins and Uestwater [54], among others, 
also experimentally showed that the product f . D^, given by Eq. (26), 
remains constant for certain liquids and for certain conditions. 

Since Eqs. (26) and (28) are cocmonly used and Indicate a reduction 
In frequency as the departure diameter Increases, both relations are 
plotted In Fig. 88, using the present experimental results obtained for 
various orientations with F113 and LN 2 . The results indicate a marked 
Increase In the f • products with increasing heat flux. This ag; , ees 


! 
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with the findings of Ref. [48]. On the other hand the f . products 
do rot vary significantly with Increasing teat flux and almost remain 
constant, which agrees with the findings of McFadden and Grassmarr [53]. 

C. EVALUATION OF NUCLEATE BOILING MODELS 

Up to this point the experimental results related to the various 
parameters such as surface temperature, active site density, average 
population density, frequency of vapor bubble emission and maximum bubble 
size have been discussed. In this section discussion Is directed towards 
theoretical models. First the separation of the natural convection and 
nucleate boiling components of heat flux Is presented. Then two theor- 
etical models, whose validity It is possible to check with experimental 
data obtained during this investigation, are considered. 

Merte and Clark [55] considered that natural convection and nucleate 
boiling act Independently of each other on a given surface, and that the 
heat transfer by nucleate boiling, Q^, and heat transfer by natural 
convection, Q N q, were confined to the related area fractions A^g/Ay, and 
A^/Ay, respectively. Therefore, the total heat transfer to the liquid 
could be written as 

t Q,mc (29) 

and the total heat flux to the liquid Qt /Ay could be written as 

(30) 

Rearranging for Q^/Qt . the nucleate boiling teat transfer fraction .gives 
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Approxlmate relations for the natural convection heat flux frac- 
tions, (Q^c/A nc )/(Q|/Aj) , are tabulated In Table (3), using appropriate 
natural heat transfer correlations for each orientation along with prop- 
erties of the two fluids used. The measured surface-fluid temperature 
differences are sifcstltuted Into the natural convection heat flux frac- 
tions given In Table (3) for the tests conducted with FIT 3 and LN^ at 
a/g=l and plotted In Fig. 89 as a function of heat flux to the liquid. 

As could be anticipated, Fig, 89 Indicates that the natural convec- 
tion heat flux fraction decreases with Increasing heat flux, since the 
active site density and frequency of bubble departure Increases with 
heat flux, thus a greater portion of the increased Imposed heat flux 
will be transfered by the boiling process. A similar trend was observed 
with a horizontal up orientation in FI 1 3 [34]. 

The natural convection area fraction A^/Ay Is obtained as In Ref. 
[34] by subtracting the total area occupied by vapor bti>bles at a given 
instant of time from unity 


A MC _ I Anb 
At ” ” At 


“ [*(-£) (d»*«.) (-J-) 


(32) 


The total surface area occtpled by vapor bubbles Is given by the mea- 
sured projected area of vapor bubbles, the average population density 
(q/A) and a constant to condensate for the area of Influence of the 
bubbles. Han and Griffith [56] reported that the area Influenced by the 
departure of a spherical bubble in the water had approximately twice the 
diameter of the sphere, which suggests a value of K»4. On the other 
hand, Ref. [34] obtains unreasonable results with K*4 and suggests using 
K»3. K-1 Is a lower limit and K*4 could be the t#per limit. In order 
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to conpare the results of the present work with th^ results of Ref. [34] 
K«*3 Is also used in the present work. 

Using K°3 in Eq. (32) along with experimental results of U an,-; 

n>aX 

n/A, A N( ,/A Is calculated. Substituting the values obtained for the nat- 
ural convection area fractions and the natural convection heat flux frac- 
tions obtained Into Eq. (31), the nucleate boiling heat flux fractions 
are determined, and plotted in Fig. 90 as a function of a dimensionless 
heat flux, defined as the ratio of the heat flux to the burnout heat 
flux for that liquid. The results in Fig. 90 Initcate that Q^/Qg in- 
creases toward unity with increasing heat flux, which is expected since 
the natural convection heat flux fraction decreases as the imposed heat 
flux increases. 

The next step Is the deterai nation of the total bubble heat trans- 

fer, Qjotai > defined as the quantity of heat transfered by the growth 
Bubble 

and departure of a single bubble. This quantity Is needed to show the 
significance of the nrlcrolay.-. contribution and to assess the adequacy 
of several nucleate boiling models. The overall nur boiling heat 
flux is obtained by multiplying Eq. (31) with Qi /Ay, „.e total heat 
flux to liquid. 

- JUS* _2i_ f 33) 

At Qg At 

Q^/At Is the total heat flux associated with the action of bubbles and 
is plotted In Fig. 91 as a function of heat flux to the liquid. It shows 
the effect of heat flux and orientation on the nucleate boiling heat flux. 
The 45° line on Fig. 91 represents the results that would have been ob- 
tained at absolute zero gravity a/g»0 condition, and Is approached quite 
closely by the horizontal down orientation, where the influence of nat- 
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ural convection Is a mini mum. Q^/Ay Increases with heat flux, and for a 
given heat flux It Is slightly higher with vertical orientations than 
with horizontal up orientations. 

Finally, dividing the nucleate boiling heat flux, Q N g/A T , by the 
active site density and the frequency of vapor bubble formation gives the 
total heat transfer associated with £ single bubble. This Is defined as 
the total bubble heat transfer and given by 


a 


toIjI = 


Q* / At 


(34) 


Figure 92 shows the variation of the heat transfer associated with 

a single bi&ble with heat flux and orientation for both F113 and LN 2 re- 

spectlvely. In contrast to the findings of [34] for saturated F113, 

Fig. 92 Indicates that Qrotal decreases with Increasing heat flux. 

Bubble 

This could be due to a more rapid Increase of N/A and f than the heat 

flux Itself. As discussed in earlier sections, both N/A and f Increase 

as power functions of the heat flux (e.g. Figs, ffj-84) while heat flux 

Itself Increases linearly. There Is no significant difference In Qrotal 

Bubble 

between the vertical and horizontal up orientations, OTotal consld- 

Bubbie 

erably higher with the horizontal down orientation. Several factors may 
contribute to this difference: \) A higher microlayer contribution with 
the horizontal down orientations, discussed previously; 2) buoyancy 
forces hold the vapor bubbles against the surface, this increasing the 
contact time of each bubble with the test surface, which results in larger 
bubble diameters and at the same time reduces the frequency of the bubble 
departure. All these effects increase the heat transfer associated with 
the Individual btAbles. 
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The results of Ref. [34] are also shown in Fig. 92 for comparison. 

The difference between the result of the present investigation and the 

results of [34] Is due mainly „o the difference in the determination of 

site densities where for a given heat flux Ref. [34] reported an active 

site density about four times as large with oxide coated glass in F113. 

The explanation for the discrepancy was given In Section B-2. Figure 92 

also indicates that a similar decrease in Qrotal with increasing heat 

Bubble 

flux occurs in LN ? as with FI 1 3 . Qiotal for LN, Is somewhat higher 

Bubble 

than with F113, which could be due to the higher latent heat of vapori- 
zation of LN 2 compared to F113. 

Now that values of Qjotjij are <teterained from, the experimental 
measurements of surface temperature, active site and average population 
densities, frequency of vapor bubble emission and maxim#! bubble size 
obtained in this present investigation. It Is possible to compare th^se 
results with the following models. 

1. Latent Heat Transfer Model 

It was first suggested by Jakob [51] that the rate of heat transfer 
in nucleate boiling Is related to the bubble voli®e of the latent heat. 
This theory has been supported by the recent experimental results of 
saturated water boiling on a horizontal wire [36]. In terms of time 
averaged measurements obtained In the present investigation the latent 
heat transfer model may be expressed as follows 


Hut 


( 35 ) 
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Q rM = = p, hf» VmHI. (36) 

BaUU ( N /AtJT ' 

The experimental results of the present Investigation and Ref. [34], 
along with results predicted by Eq. '36) are presented In Fig. 93 for 
comparison. According to Fig. 93 the latent heat transfer model predicts 
values considerably lower, Indicating that other Mechanisms of heat trans- 
fer are acting, such as pumping action and Increased forced convection 
due to change In the momentum of liquid panicles as a result of bubble 
growth and departure According to Fig. 93 the horizontal down orienta- 
tion gives results closest to those predicted by Eq. (36). With the 
horizontal down orientation the superheated liquid layers remain near the 
surface, and the pumping and forced convection actions by the bubbles be- 
come less effective. Therefore, the contribution of the latent heat 
transport to thet total Is larger with the horizontal down orientation 
than with the vertical and horizontal up orientations. 

2. Vapor Liquid Exchange Model 

In this model, postulated by Forster and Grief [37], heat transfer 
by nucleate boiling was described as the exchange of a bifcble volune of 
liquid at tenperature Ty for a bubble volume of liquid at the bulk tem- 
perature Ti due to the action of bubble departure and collapse expressed 
by 

-*)(.£ )f 


(37) 
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or 

a ISli'l. /vs«kW« = $^7%~ = (1 C « (Tt - T 0 (38) 

Figure 94 presents the experimental results of Qjotal along with 

Bubble 

the theoretical predictions computed from Eq. (38). The results of Ref. 
[34] are also Included for comparison. 

As Fig. 54 indicates, the vapor liquid exchange model largely over- 
predicts the heat transfer per bifcble. A similar conclusion Is reached 
in Ref. [34]. Part of the deviation could come from the fact that the 
liquid near the surface is not uniformly at the wall temperature, as 
the model assumes. The deviation is largest with the horizontal down 
orientation and smallest with the vertical orientations. Due to less 
mixing action the temperature gradient may be expected to be hicfrest 
within the liquid with the horizontal down orientation, resulting in the 
large deviation observed. 



CHAPTER VII 


CONCLUSIONS 

1. This research has demonstrated that semi-transparent thin gold 
films, properly calibrated, can be used simultaneously as a heater and re- 
sistance thermometer for temperatures ranging from -320°F to +25Q°F, with 
better than -2°F accuracy and producing heat flux levels up to burnout in 
LN^ and F113. The resistance of the gold film surfaces increases with 
time, but the slope dR/dT remains constant. Therefore, once a complete 
calibration of a test surface has been made. It is sufficient to make oily 
a single point calibration prior to each test run. 

2. From the step increase of power input up to the initiation of 
boiling two distinct regimes of heat transfer mechanisms exist; the con- 
duction and convection dominated regimes. The time duration in each regime 
is considerably shorter with LN 2 than with F113, and decrease as heat flux 
increases, as gravity Is reduced, and as the orientation is changed from 
horizontal-up to horizontal -down. These changes are a consequence of 
changes In boundary layer temperature levels associated with changes in 
natural convection. 

3. With a step increase of power input, boiling initiated first at 
a single point and then spreads over the surface. The delay time for the 
inception of boiling at the first site decreases as heat flux increases, as 
gravity is reduced and as the orientation is changed from horizontal up to 
horizontal down. The delay time depends on the transient growth of a 
thermal boundary layer. Therefore, as heat flux is Increased, the liquid 
film tenperature increases more rapidly. Similarly, with a reduction of 
gravity or changing orientation from horizontal up to horizontal down, con- 
vection effects are reduced or eliminated, and more energy is available 
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to heat the boundary layer. The delay time also defends on the liquid 
properties, being shorter with LN^ than F113. The spread velocity of 
nucleate boiling increases with increasing heat flux, and with a reduction 
of gravity. It is larger with LN 2 than with F113, and is larger with 
smoother surfaces. 

4. With a stepwise increase in power input, a heater surface tem- 
perature overshoot and undershoot occur before the steady state boiling 
temperature level is established. The degree of overshoot varies depending 
upon the heat flux, gravity, orientation and surface characteristics. It 
is higher with smoother surfaces and increases with increasing heat flux. 

The effect of orientation on the degree of overshoot depends on gravity 
and the heat flux level. At lower heat flux levels the highest overshoot 
occurs with vertical surfaces, while at the higher heat flux levels, it is 
highest with the horizontal up orientations. It is lowest for horizontal 
down orientations for all heat flux levels. 

5. The surface superheat at the inception of boiling is Independent of 
heat flux and orientation. It is a function of gravity, however, being 
smalle* at a/g=o than at a/g»l. It also appears to be a function of sur- 
face roughness. 

6. While the active site and average population densities increase as 
a power function of the heat flux, these are Influenced by orientation and 
gravity to a minor extent, only. 

7. Experimental results indicate that the maximum bubble size and the 

frequency of bubble departure are correlated for various levels of heat 
flux and orientation by the relation, (D f = constant. 

3. The heat transfer associated with each departing vapor bubble, 

QTotal , decreases as heat flux Increases. 

Bubble 
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Fig. 4. Accelerometer measurements with free fall 
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Fig. 9. Thin film surface dimensions 
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Fig. 12. Calibration test, gold film PCG 6 
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Fig. 13. Calibration and power circuit 










Fig. 14. Instrumentation and power supply 
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Fig. 15, Sequencing and delay circuit 
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Fig. 16. Heat loss through glass substrate, FI 1 3 
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Fig. 17. Heat loss through glass substrate, IN 
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Flo. 29. Effect of heat flux, orientation and gravity on ten^erature overshoot 
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Fig. 33. Comparison of computed and measured transient 
thin film temperature, LN ? 







Fiq. 35. Onset of convection, LN 9 , horizontal -up, a/q 
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Fig. 41. Elapsed time for Incipient boiling, and maximum surface 
temperature, F113, verticil, a/g-1 
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Fig. 46. Elapsed time for spreading of boiling, T« F113, 
horizontal-up, a/g»l and « 0 5t 
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Fig. 56. Elapsed tine to establish steady state nucleate 
boiling FI 1 3 and hori zontal-down , a/g=l 



Fig. 57. Elapsed time to establish steady state nucleate 
boiling F 11 3 and LN ? , all orientations, a/g - 0 
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Fig. 92 . Total bubble heat transfer 
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APPENDIX A 


ONSET OF CONVECTION 

With the transient process being used here, a step Increase In hLat 
generation In the thin film deposited on t\a glass surface Is used In order 
to achieve a step Increase In heat flux to the liquid. The quantity mea- 
sured Is the heat generation In the film, which Is distributed to Doth 
the liquid and the solid. The fraction going to the liquid can be com- 
puted If it Is possible to calculate the fraction going to the solid siae. 

A solution to this problem Is available provided that the glass plate can 
be considered as a semi-infinite solid during tiie time period of Mterest. 
To assess this, use Is made of the expression for the penetration depth 
of temperature with an Imposed heat flux [57] 

c5 = (6at )'* (39) 

Substituting the properties of pyrex Into Eq. (39), it will take 
approximately 6.0 seconds for the far side of a 3/16 inch thick plecs 
of pyrex to start rising In temperature. The maximum time periods being 
considered In our experiments are approximately 2.5 seconds, so the py- 
rex plate can be considered as a seml-lnflilte solid during this period. 

The problem whose solution Is now desired Is that of two semi- 
infinite solids with a step Increase In a plane heat source at the boun- 
dary of these two solids. The pyrex substrate Is one solid, and the 
liquid being heated represents the other solid, up to the point where 
convection or boiling begin. The heat capacity of the thin film can be 
neglected. The solution Is presented In Ref. [58], and reproduced below: 
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T / , T. _ 2 (fiA)r('i^i) ’ r 


(40) 




X. 

2(**l Y 


(41) 


For a plane heat source of constant strength Q/A , It can be shown, 
with Eqs. (40) and (41), that the ratio of the heat flux to each region 
Is constant and given by: 


^ ) /l (42) 

<3, k, \ di J 

The fraction of the total input which goes Into the liquid Is then given 
by 


[ a / A Jf , ! (43) 

( a AX i + V<s t 

The temperature distribution In the liquid (x,t) is of interest 
In establishing the conditions for Incipient boiling under the transient 
process. The Interface temperature, Ty , corresponding to X»0 In Eqs. 
(40) and (41), Is 


T T 2 ( Q /A )t (^3 Y* ft 
T ~ 

Equation (44) can also be rewritten as 


(44) 


Tr-Tj 

(3/a)* 


_ 2 

-It* 



( 45 ) 
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Subs ti tuti ng the properties of liquids and pyrex at the corresponding. 

temperatures into Eq. (43), the fraction of the total energy input 

transferred to the liquids are: 

To FI 13 _ 


(Q/a )j 

( 6LM ) T 


- = 0 166 


— 1 f 1 1 s 


(46) 


To LN 2 


(S/aU 


|= 0-491 


(47) 


w 

Substituting properties into Eq. (4S) and rearranging will give the 
temperature rite of the gold film: 


In FI 1 3 


In LN 2 


Tr ~ Tj 

(a/A ) 4 


= 0.0223 \! 1 

FII3 


(48) 


= 0 - 00666 L ( 49 ) 

The next process which occurs in a gravitational field prior to the 
inception of boiling Is the onset of natural convection. Onset of nat- 
ural convection has been characterized as an instability problem, and 
has been reasonably well described by Vest and Lawson [T59] for the case 
of a horizontal wire. 

Rayleigh [150] showed that the stability of a horizontal layer of 
fluid heated from below Is governed by the numerical value of the d1»- 
ensionless parameter 


Tt-T; 

(a /a)i 



- 


(50) 
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It has been shown by Jeffreys [61 J, Low [62], Pellew and Southwell [ 63 ] 
that the critical Rayleigh number Is about 1700 for a fluid between two 
solid boundaries. Low [62], Pelle* and Southwell [63] calculated the 
critical Rayleigh nunfcer for rigid lower and free upper boundaries, and 
for constant wall temperature. They found that the critical Rayleigh 
number, for such a fluid layer with constant underside wall temperature 
is about 1100. Sparrow et al [64] showed that the critical Ra*800 for 
fluid layers heated from below with constant imposed heat flux. From 
Eq. (50), it is then possible to compute the thickness which becomes 
unstable. & can be expressed as a function of time using the approxi- 
mate solution of the conduction problem In a semi-infinite solid with 
an imposed constant heat flux, expressed In terms of a "penetration 
depth" gieen by Eq. (39). The approximate solution also gives 


T t (±)- T , = (_i_j 


['*■ ( a/A h (a-l-f 2 


kl 


(51) 


The temperature gradient within the penetration depth, assuming it to 
be linear, is given by 


yg __T vlihJL - 

' ~ 5 


3TT 


\ o ./ a \ 

Hi 


(52) 


Substituting Eqs. (39 ) and (52) into Eq. ( 50 ) , and taking the critical 
Ra number as 800, one obtains 

(<Va)i Ys a 1 

kx7 




r, - 7 


(53) 


where 7 ] Is the time at which the instability sets in as manifested by 
the onset of natural convection. 
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Substl tutlng the properties for LN 2 Into Eq. (53), this reduces to: 

T, = 13-1 (54) 


( a/A \ 
For FI 1 3 It reduces to 




r, _ is 

'fllJ 

Table 5 gives the simary of the results similar to Eq. (54) ob 


(55) 


tained for LN^ and F113 using the three different critical Rayleigh num- 
ber's. 


TABLE 4 


The Critical Raorlel gh Number and Time 
Delay for Onset of Convection 


Liquid 

Type of 

Boundary Condition 

Critical R a 

[(Q/A) ] ,/2 r, 

ln 2 

Solid-solid 
(constant wall temp) 

1700 

19.08 

FI 13 

Solid-solid 
(constant wall temp) 

1700 

37.5 

ln 2 

Sol id- free 
(constant wall te^>) 

1100 

V>.35 

FI 13 

Solid-free 
(constant wall temp) 

1100 

24.4 

ln 2 

Solid-free 
(constant heat flux) 

800 

13.09 

FI 1 3 

Solid- free 
(constant heat flux) 

800 

18 


The above development applies to the onset of natural convection on 
a horizontal heated surface facing upward. The onset of natural convec- 
tion on a vertical heated surface is treated some^at differently. 
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As first pointed out by Slefel [68]. the Initial behavior of the 
temperature and velocity field for a semi-infinite vertical flat plate 
are the sane as for the doubly Infinite vertical plate. The temperature 
field could be obtained by solving an unsteady one dimensional heat- 
conduction problem. The transition to convection begins only when some 
effect from the leading edge has propagated up the plate to the parti- 
cular point In question. Before this leading edge effect has propagated 
to a particular point the region effectively does not know that the 
plate has a leading edge. Goldstein and Briggs [65] obtained an equa- 
tion for predicting the penetration distances of disturbances due to the 
Influence of the leading edge for a vertical semi -Infinite flat plate, 
for boundary conditions of a step change In surface heat flux, based on 
the velocity and temperature fields along a doubly Infinite flat plate. 
The maximum value of the penetration distance, Xp, obtained at any time 
for arbitrary Prandtl number was 

X _ •oossir(*/A)i{«*-) /k t l (56) 

Substituting the height of the test surface for Xp and liquid prop- 
erVias, the following relations are obtained: 

For FI 13 

( a/4 iL Ti = 21 - 3 <57) 

For 

( r , = 16.1 <*» 


a 
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values In Eqs. (57) and (58) give the time required for the 
leading edge effect to reach the apper edge of the vertical test sur- 
face, which represents in turn the completion of the onset of free con- 
vection. 



APPENDIX B 


ERROR ANALYSIS 


In the following section the errors associated with the 
determination of test surface temperature, area, and heat flux are dis- 
cussed. To obtain the percent of uncertainty in the quantities calcu- 
lated the procedure described by Kline and McClintoc [66] Is used. The 
uncertainty associated with computation of the result R, which is ob- 
tained by combining certain numbers of variables , V 2 , Vg...V n , 
given by 


W/9 


(BR Wl f , 

(a* W a\ z . 

Id R 

Wn f 

R 

r*. 

\[ *7, -*l + 

(all R ) + ' 

"'VWn 

R / 


where Is the uncertainty In the result R and W^. ..W n are the 

uncertainties associated with each Independent variable. 


A - TEST SURFACE TEMPERATURE 

To calculate the uncertainty associated with the determination of 
test surface temperature, the uncertainties associated with Ry, R , T , 
dT/dR must be determined, as is noted in Eq. (5y). 


1) Test Surface Resistances 

To oteasure the test surface resistance the set up shown In Fig. (13) 
is used. Ry is given by 


D — Er 
Kt ” Es 


The uncertainty associated with Ry Is given by 

r- - * •n'fc 


__ 
Rt 


(t T f + 


(60) 


(61) 
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For calibration runs, Ey and E $ are measured wtth a K-3 potentiometer. 
R $ is the .50729 31 standard resistance. The uncertainties associated 
with those are W^y/E . 0 1 5t , W^; s /Ey« .015% and WR<./R s “.Q25t. Combining 
all these we obtain the uncertainty in computing R c to be 



For regular runs, Ey and E $ are obtained from Sanborn recordings. 

The uncertainties due to non-linearity of *.25 division in the recorder 
amplifier are W£y=>12 .5x10'^ volt and W£«j*6.25x10" 6 volt at corresponding 
sensi tl ti ties of .005 volts/div. and 50 mv/div., respectively. The un- 
certainty associated with the .01J1 standard resistance is *.04?. Fora 
typical case Ey=10.8275 volt and E s -20.815mv related to percent uncertain- 


ties were, W^/E s =.03% and W£y/E-p“ .01 16X. Substituting these values in 
Eq. (62) gives , 

- 1 - = (M6)%4%3 2 v| 6*=. Q51% (63) 

J 


2j Calibration Temperature T c 

To obtain calibration temperature the E.M.F. generated by a Cu- 
constantan thermocouple was measured by the K-3 potentiometer. The un- 
certainty associated with the K-3 potentiometer is +0.15% + G-*/ xV 


For B. c = /970 jm\t W Ec » i. 1Z8 


(64) 


and using Leeds and Northrup standards, the percent incertainty Is cali- 
bration temperature, T c , was found to be Wj c /T c “.008? for Freon 113. 
Following similar procedures and using NSD standards, the uncertainty 
in T c for LN 2 was found to be *.036t. of 139°R which is less than *.05°R. 
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3) Calibration Curve Slope dR/dT 

In determination of calibration curves for the test surfaces at FI 13 
saturation temperature levels, a Mueller Brldgeand platinum resistance 
thermometer were used. The uncertainties in the Mueller Bridge readings 
were *.021. The uncertainties in the platinum resistance thermometer 
arises from two sources: 1) the error In the variation of R Q which is 

*.02°C, 11) uncertainty coming from the assumption of linear variation 

between tabulated values which is *.001°C. The first error was always 
Included In calculations. Therefore, during the calibration of a test 
surface, surface resistances were measured with an uncertainty of 
W RW /RW n i.0002 and liquid temperatures were determined with an uncertainty 
of *.01°C. Thus the two ends of the calibration curve have the follow- 
ing values 

T, = 117. 5 + 0016’F 65 .25*.00(8> 

(65) 

R,= 5.0844- 001 R 2 = 4.949+ 00099 

For the worst case the uncertainty In the slope dT/dR is *1.46%. 

For calibration of the test surfaces in LN^, the set up shown In 
Fig. (1) of Ref. [41] Is used, along with copper- cons tan tan thermocouples. 
Test surface resistance is obtained with an uncertainty of *.033% and 
liquid temperature Is determined with an uncertilnty of *.036%. 

Following the same procedure used for F113, calibration curves of 
the test surfaces In LN 2 are determined with an uncertainty of *2.6% 

In their slopes. 

Test surface wall temperatures are calculated from Eq. (11) 

Tf = ■£ +• m (R T - R c ) 


(ID 
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Tne uncertainty associated with Ty is obtained from the following 
equation 


Wr T 


f z 2 4 V* 

' (Wr t 'f + + (W*-«-c) j> 


(6b) 


combining all these, gives the uncertainty in computing wall temperature 
of the test surface in FI 1 3 to be -1.33°F. Following the same proced- 
ures and using corresponding uncertainties, the test surface temperatures 
in boiling with LP^ are determined *1.68°R uncertainty. 

B. AREA MEASUREMENTS 


Uncertainties associated with the determination of the portion of 
the test surface area photographed arises from 1) uncertainty in mea- 
surements of the actual distance between dimension lines, 2) uncertainty 
in measuring actual test surface area, 3) uncertainties In measuring 
linear distances on the analyzer screen, 4) uncertainty in tracing the 
area from the analyzer screen, 5) uncertainties In determining enlarged 
areas by means of a planlmeter. 

The spacing between the dimensioning wires and the test surface 
dimensions are measured by a .001" traveling microscope. Therefore, un- 
certainty associated with .1 spacing and test surface area are *.15 and 
.14% respectively. Linear dimensions and the areas measured by hairlines 
on the analyzer screen have an uncertainty of -.55 and *.75 respectively. 
Tracing enlarged areas from the analyier screen Introduces an uncertainty 
of *.55. Finally determining traced areas by using a planimeter intro- 
duces an uncertainty of -.65. Therefore, the uncertainty in determina- 
tion of the areas Is 

’ / a 

(£<♦)* + (7of+(50) l 4(60) Z M0*=+1.06/J (67 ^ 


- + 



! 

I 
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C. HEAT FLUX 

The uncertainties In determination of heat flux arise form a) the 
uncertainty associated with the determination of total heat flux from 
the relation 


(«/a) 




EtI 

A 



( 68 ) 


and b) the uncertainty associated with heat loss from the test section. 


a) The uncertainties associated with the determination of the 
total heat flux Is given by 


... 

(Q/a) 




(69) 


where 


W*r 

Er 


im%i = = &- = *l.0bZ (70) 

bs A 


Substituting these values In Eq. (69) will give 


W(< a A) . 

{G/a) ~~ 


(71) 


b) Uncertainties associated with the heat loss from the test sec- 
tion Is directly related to the uncertainties associated with the deter- 
mination of (T-p-T^) ( 


1A/(T t -T|) 


(w Tr f + (Wt<)‘ 


(72) 


where W^»1.337, W^«.061 as found before. Substituting In £q. (72) 
gives W^_y )»1.338®F which In turn gives *2.1% uncertainty In heat loss. 



APPENDIX C 


PHOTOGRAPHIC SAMPLES 

In this section typical samples of photographic and thermal data 
are presented along with wie results obtained from the analysis of these 
data. 

Figures 95 and 96 Illustrate typical photographic and thermal re- 
cordings obtained from transient boiling on a vertical surface In F113 
at a/g-1 and a/g»0, respectively. These figures demonstrate the Incep- 
tion and spreading of boiling on a vertical surface at standard and zero 
gravity. The corresponding variation of the test surface temperatures 
are also plotted In these figures. 

Figures 97 and 98 are typical photographic samples of steady state 
nucleate boiling for a vertical surface In FI 1 3 at a/g»l and a/gO, re- 
spectively, and were obtained consecutively by starting the camera well 
before the test package was released. The pictures were all taken at a 
framing rate of 2200 pictures per second. Figure 97 shows the sequence 
at a/g*l beginning 0.13 seconds prior to release, and covers the time 
period of 0.02 seconds. All photographs obtained In this Interval are 
reproduced In Fig. 97. Figure 98 shows the sequence at a/g*0, beginning 
0.15 seconds after release of the test packige. The samples reproduced 
In the upper row are for each tenth frame obtained In the time Interval 
of 0.15 seconds to 0.20 seconds after release of the test package, 
while the samples In the lower row are for each tenth frame obtained In 
the time Interval of 0.25 to 0.30 seconds after release. 

The high speed pictures of boiling clearly show that bubbles con- 
tinue to grow and depart from the active site during the early part of 
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the 0-G. However, during the later part of 0-G they remain stationary 
on the surface and grow without departing. Figures 99-101 Illustrate 
this phenomenon. In Figs. 99 and 100, data for the same site Is given, 
but with different time scales. Figure 99 shows the growth of a bubble 
at about .02 seconds after the Inlt.atlon of 0-G. The diameter of the 
bubble was about 0.001 Inches when first observed. It then grew to a 
diameter of about .03 Inches In .04 sec., then departed from the surface. 
After about .15 seconds elapsed In 0-G, another bubble appeared at the 
same site but did not depart from the surface, but rather continued to 
grow until It coalesced with ether bubbles. Thl oble appeared at the 
site .146 seconds after the release of the test * -r-.'.qe, with a diameter 
of .0013 Inches and grew to a diameter of .084 inches until it coalesced 
with other bubbles growing at neighboring sites at about .34 sec. after 
test package release. Figure 100 shows the sequential growth of these 
two bubbles on a common time scale. The time scale of Fig. 100 Is 50 
times larger than that of Fig. 99. 

The departure of the bubble In the early stages of 0-G Is believed 
to be caused by the Inertial effects of the liquid motion set up during 
1-G boiling. After a certain time in 0-G the liquid motion dies out due 
to viscous effects and the bubbles remain on the surface. 

Figure 101 shows the bubble growth at 0-G for various heat Fluxes 
after the liquid motion has died out. 
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Fig. 97. Typical steady state nucleate boiling, HI* vertical 
a/g=l, run 44, Q/A=4200 Btu/hr 
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TIME -FRAMES (at 1080 Frames / Second ) 

Fig. 101. Typical bubble growth for various heat fluxes, a/g 



APPENDIX D 


CALIBRATION RESULTS OF THIN GOLD FILMS 

In this section calibration curves of various pyrex test surfaces 
coated with thin gold films are presented to show the type of shifting 
occuring In the resistivity of gold films with time and the effect of 
tempering on shifting. 

Figures 102-105 are the calibration curves of the test surfaces 
PCG-15, 16, 19 and 20. Figures 1% and 107 are the calibration curves 
of PCG-21. Figures 108-110 are the calibration curves of PCG 22-24. 
Figures 111 and 112 are the calibration curves of PCG 29. Figures 102, 
104-108 and 112 are the results of calibrations at LN 2 temperature level. 
Figures 109-111 are obtained at F113 temperature level. Figure 103 
covers the complete temperature range. 

Test nifflbers on the calibration curves Indicate the sequence of 
each data point determined during calibration. Each data point repre- 
sents mean of four measurements with currents between 5-30 oa . which 
corresponds to test surface heat fluxes, depending on the resistivity 
of the test surface of between .01 to 1.0 Btu/hr-ft^. Calibration dates 
are also given on the figures. 

The results presented In Figs. 102-112 clearly show that the slope 
dR/dT always remains constant Indicating that once the complete cali- 
bration of a test surface had been made It was deemed sufficient to 
make only a single point calibration prior to and following each test 
run. 

All calibration curves presented in Figs, 102-112 except Fig. IDS 
are obtained after the test surfaces are baked for about 16-18 minutes 
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at approximately 275°C. Figure 106 Is obtained before tempering. Com- 
parison of these two shows the effect of tempering on thin gold films. 
In this case the resistivity of the test surface PCG-21 at 138°R was 
dropped sharply from 1.055 to ,88 A/SQ. 

Test surfaces PCG 15, 16, 19 and 20 were kept In nitrogen with 3i 
H 2 at one atmosphere during the baking, and surfaces 21-24 and 29 were 
tempered in atmospheric air. No effect of oxidation during tempering 
is noticed. 
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APPENDIX E 


TEST SURFACE TEMPERATURE 

The test surface temperature as measured Is a mean value, and is 
determined from the total resistance of the surface. The surface resis- 
tance at any tine is equal to the slot of the resistance of boiling part 
(Rq) and the resistance of the nonboiling part (R^g) of the surface. 



Assuming that boiling spreads linearly across the surface as shown 
above, the total surface resistance can be approximated as 

Rt — Rg •+ Rwg (73) 

Kor a given width and surface thickness, the resistance can be written 
as a function of the length of the surface as 

R r 2 p a XU)+£, [l-X({)] (74) 

Since the average surface te^>erature is a linear function of total 
resistance 
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T » R( 


(75) 


Equation (27) can ilso be written as 


Rf — Cont4jnf T 


(76) 


where the constants in Eqs. (75) and (76) are different. Equation (76) 
■ay be written as follows 


R t = Consijnl J TaXU) + Twe(f) L~X(t) 


(77) 


Since Tg is not a function of time differentiating and equating the 
Eq. (77) to zero will give the condition for maximum surface temperature. 

[r„ui-T»] =. [l-xco]^*! m 

Multiplying both sides of Eq. (78) by dt and separating the variables 
ATu*(i) dKU) 


T*a(0-Ta l-Ui) 

Integrating Eq. (79) will give 


(79) 


jjVi»(0- Taj Ij.-XU'jJ* Convbnl 


(80) 


where X(t) 1 L and (t) f T g 

Multiplying both sides of Eq. (82) by the width of the test sur- 
face will result following relation 


pTwaU)- T»J [a t " A*J ■» 


(81) 


Equation (81) is an approximate relation and the constant In the 
equation need not be evaluated. This relation can only be used to support 



qualitatively, the experimental results obtained In present Investiga- 
tion which showed that surface temperature continued to Increase even 
after inception of first site, until a certain percent of the test sur- 
face area Is covered by the active nucleating sites. Equation (81) Is 
hyperbolic function, which indicates that as Ay-A B decreases T NB -T B 
Increases. 
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